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The  turbulence  characteristics  of  a  water  jet  flowing  near  a 
free  surface  were  studied.  The  focus  of  the  project  was  to  develop 
two-point  closure  methods  which  may  make  possible  a  "universal" 
description  of  the  turbulence  in  the  region  where  the  jet  interacts 
with  the  free  surface. 

The  interactions  of  a  planar  jet  with  an  air-water  interface 
were  determined  by  measuring  turbulent  velocity  fluctuations  and 
pressure  fluctuations  using  hot-film  anemometry  techniques ,  A  two- 
dimensional  jet  was  formed  by  pumping  water  through  a  rectangular 
slit  into  a  channel  filled  with  still  water.  The  raw  data  are 
sampled  time  records  of  the  anemometer  output  voltages  obtained  at 
different  positions  in  the  jet.  These  time  records  together  with 
the  sensor  calibration  information  are  analyzed  using  digital  signal 
processing  techniques.  The  double  and  triple  velocity  correlation 
functions  and  pressure-velocity  correlation  functions  as  well  as  the 
mean  velocities,  Reynolds  stress,  and  turbulence  intensities  were 
calculated. 

These  calculated  results  were  used  to  develop  a  two-point 
turbulence  closure  scheme  for  the  jet  flow.  The  closure  model 
contains  the  effects  of  small  non-homogeneities  and  anisotropy, 
since  there  is  a  preferred  direction  in  the  jet  flow.  Because  of 
the  proximity  to  the  free  surface,  the  flow  field  is  not  symmetric 
about  the  axis  of  the  slit.  By  means  of  a  projective 
transformation,  the  jet  flow  turbulence  results  can  be  used  to 
predict  the  turbulence  in  the  zone  of  interaction  between  a  wake  and 
a  free  surface. 

SUMMARY  OF  THE  EXPERIMENTAL  WORK 


The  chronology  of  the  project  is  as  follows.  During  the  first 
year,  the  flow  system  was  designed  and  built,  velocity  measurement 
techniques  were  developed,  and  velocity  sensor  calibration  equipment 
and  methods  were  developed.  The  unique  feature  of  this  calibration 
method  was  devising  a  systematic  method  for  dealing  with  sensor 
fouling  during  long  periods  of  continuous  operation  of  the  hot-film 
sensors.  This  work  was  the  subject  of  the  M.  S.  thesis  written  by 
Eric  J.  nine.  During  the  second  year,  techniques  for  measuring 
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turbulent  pressure  fluctuations  were  devised.  A  commercially 
available  sensor  was  used,  but  extensive  development  work  was 
required  to  adapt  it  to  operation  in  water.  This  sensor  system  was 
used  along  with  a  standard  hot-film  anemometer  system  to  measure 
two-point  pressure-velocity  correlations.  These  results  were  the 
subject  of  the  M.  S.  thesis  written  by  Lisa  R.  Harris.  During  the 
third  year,  the  techniques  developed  were  applied  to  an  extensive 
study  of  the  jet  flow.  The  mean  velocities  and  jet  boundaries  were 
measured  along  with  single-point  correlations — turbulent  intensities 
and  Reynolds  stresses.  Two-point  correlations  were  also  measured. 
The  single-point  and  two-point  measurements  are  the  subject  of  the 
M.  S.  thesis  written  by  Eric  Menning. 


SIGNIFICANT  EXPERIMENTAL  RESULTS 
A.  Development  of  Experimental  Techniques 

Sensor  Calibration  One  important  accomplishment  has  been  to 
develop  a  method  for  allowing  for  fouling  of  hot-film  sensors  when 
making  measurements  in  untreated  water.  Fouling  cannot  be  avoided, 
but  it  has  been  shown  that  the  decrease  of  anemometer  output  voltage 
with  time  while  making  measurements  in  fouling  conditions  is  linear 
except  for  an  initial  "seasoning"  period.  By  calibrating  the  sensor 
periodically  during  the  course  of  a  series  of  experiments,  and  by 
keeping  careful  account  of  the  time  at  which  individual  measurements 
were  made,  a  sensor  calibration  relationship  for  each  individual 
measurement  can  be  calculated.  The  sensor  calibration  was  done 
using  a  laminar  jet  apparatus  supplied  with  the  same  untreated  water 
which  was  being  used  for  the  jet  flow  experiments. 


Pressure  Fluctuation  Measurement  Another  important  accomplishment 
has  been  the  successful  application  of  a  bleed  type  pressure  sensor 
to  measuring  turbulent  pressure  fluctuations  in  water.  Until 
recently,  pressure  sensors  sensitive  enough  to  measure  turbulent 
pressure  fluctuations  were  not  available,  so  the  effect  of  the 
pressure  terms  could  not  be  included  in  two-point  turbulence  closure 
models  for  nonisotropic  turbulence.  Jones  (1981)  and  Spencer { 1970 ) 
report  the  development  of  a  bleed-type  pressure  sensor  based  on  a 
hot-film  anemometer  system.  A  hot-film  velocity  probe  connected  to 
a  constant  temperature  anemometer  circuit  is  enclosed  in  a  tube  open 
to  the  pressure  to  be  measured  at  the  downstream  end.  The  upstream 
end  is  connected  to  a  constant  pressure  source  of  fluid,  and  there  _ 
is  a  constant  flow  of  fluid  from  the  constant  pressure  supply,  past 
the  velocity  sensor,  and  into  the  flow  field  where  the  pressure  is 
being  measured.  The  velocity  of  the  fluid  flowing  through  the  tube 
is  related  to  the  pressure  difference,  so  the  downstream  pressure 
can  be  related  directly  to  the  anemometer  output.  It  is  well-known 
that  hot-film  velocity  probes  are  very  sensitive  to  changes  in  the  - 
cooling  velocity  in  the  neighborhood  of  the  sensor.  This  means  that 
pressure  measurements  based  on  the  velocity  determined  using  a  hot- 
film  sensor  are  equally  as  sensitive.  0/3 
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Pressure  sensors  such  as  described  above  are  available 
commercially  from  TSI,  Inc.  These  pressure  sensors  were  originally 
developed  for  pressure  measurements  in  gas  flows  using  helium  or 
nitrogen  as  the  bleed  fluid.  Such  sensors  had  never  been  applied  to 
pressure  measurements  in  liquids,  so  some  development  work  was 
necessary.  The  bleed  fluid  must  be  non-conducting  to  prevent 
corrosion  of  the  sensor,  and  the  flow  through  the  sensor  must  be 
laminar.  For  measurements  in  water,  the  bleed  fluid  selected  is 
ethyl  alcohol  which  is  miscible  with  the  main  water  flow.  The 
ethanol  must  be  supplied  at  constant  pressure  to  the  pressure 
sensor.  In  the  system  devised  for  doing  this,  ethanol  flows  to  the 
sensor  from  a  reservoir  which  is  pressurized  with  nitrogen.  The 
ethanol  is  pumped  by  a  gear  pump  into  the  reservoir  from  a  small 
supply  tank.  There  is  a  filter  in  the  line  to  remove  any  particles 
which  may  be  in  the  ethanol.  The  pressure  is  maintained  constant  by 
a  system  of  bleeds  and  pressure  regulating  valves. 

To  calibrate  the  pressure  sensor,  the  sensor  is  placed  in  still 
water,  and  the  bridge  output  is  recorded  as  the  pressure  upstream  in 
the  reservoir  is  changed.  Since  the  flow  rate  through  the  sensor 
depends  on  the  pressure  difference,  it  does  not  matter  whether  the 
upstream  or  downstrecua  pressure  is  changed  during  calibration. 

The  pressure  fluctuation  measurements  required  an  exacting  and 
delicate  experimental  technique,  and  the  sensor  was  rather  fragile. 
The  sensor  failed  twice  during  the  course  of  the  work  and  had  to  be 
returned  to  the  manufacturer  for  repair.  In  November  1987,  it 
failed  again  and  could  not  be  repaired.  This  means  that  we  never 
were  able  to  obtain  data  suitable  for  some  of  the  calculations  we 
wanted  to  do. 

B.  Single-Point  Results 

The  single-point  results  were  obtained  using  a  crossed-sensor 
probe  so  that  both  longitudinal  and  transverse  fluctuations  could  be 
calculated.  The  results  presented  were  obtained  by  calculating  the 
mean  of  time  series  records  which  had  been  converted  to 
instantaneous  velocity  values  by  applying  the  appropriate 
calibration  equation.  These  single-point  quantities  are  presented 
as  functions  of  position  in  the  jet.  The  results  presented  are  not 
synoptic — measurements  were  ruade  at  different  times  in  a  stationary 
turbulent  field.  The  complete  results  are  included  in  the  thesis 
written  by  Eric  Menning  which  is  attached  as  an  appendix  to  this 
report.  The  highlights  of  the  work  are  summarized  below.  The 
notation  used  to  describe  the  system  is  shown  in  Figure  1. 

Mean  Velocities  The  mean  velocity  profiles  at  different  positions 
downstream  from  the  source  for  a  submergence  ratio  of  H/D  =  6.9  are 
shown  in  Figure  2.  The  dimensional  data  show  that  the  jet  decays 
rapidly,  and  distortion  due  to  the  presence  of  the  free  surface  is 
evident.  The  center  line  of  the  jet  is  displaced  toward  the  surface 
as  the  distance  downstream  increases.  At  the  point  farthest 
downstream,  the  profile  is  not  well-defined.  The  jet  boundary  seems 
to  have  intersected  the  surface  at  approximately 


4 


x/D  =  25.  The  mean  velocity  data  are  summarized  in  dimensionless 
form  in  Figure  3  where  the  mean  velocities  have  been  made 
dimensionless  by  dividing  by  the  maximum  velocity,  and  the  vertical 
position  has  been  made  dimensionless  by  dividing  by  the  length  scale 
defined  in  the  next  paragraph.  This  figure  shows  that  the  jet 
decays  in  a  self-similar  way,  because  the  profiles  nearly  coincide 
for  all  axial  positions  shown.  The  distortion  due  to  the  presence 
of  the  free  surface,  and  the  displacement  of  the  jet  center  line 
toward  the  free  surface  are  evident  in  this  figure.  This 
displacement  is  plotted  in  Figure  4  which  shows  that  in  terms  of 
dimensionless  coordinates,  the  displacement  of  the  jet  center  line 
approaches  an  asymptote  as  x/D  becomes  lai*ger.  The  effect  of  the 
proximity  to  the  surface  on  the  mean  velocity  profile  at  a  single 
axial  position(x/D  =  11.0)  is  shown  in  Figure  5.  The  jet  centerline 
is  displaced  toward  the  free  surface  from  the  nozzle  centerline,  and 
the  submergence  of  the  jet  seems  to  have  a  small  influence  on  this 
displacement.  Velocity  measurements  close  to  the  surface  are  rather 
uncertain  because  of  waves  present  as  the  jet  intersects  the  free 
surface. 

The  decay  of  the  centerline  mean  velocities  in  a  turbulent  jet 
is  described  by — 

(Uj^/Uq)  (x/D)0*5  =  (2/4r 

c  is  Reichardt's  constant  which  must  be  determined  experimentally. 
Mih(1975)  reports  c  =  0.140.  Equation  (1)  is  plotted  in  Figure  6 
along  with  the  measured  centerline  mean  velocities  obtained  in 
several  different  jet  flow  experiments  at  different  submergence 
ratios.  The  velocities  plotted  are  the  velocities  at  the  nozzle 
centerline.  Because  of  the  distortion  due  to  the  presence  of  the 
free  surface,  the  jet  centerline  may  not  coincide  with  the  nozzle 
centerline.  This  effect  could  cause  some  of  the  data  scatter 
observed.  There  is  general  agreement  with  the  work  of  Mih  as 
represented  by  Equation  ( 1 ) . 

Jet  Spreading  The  distance  from  the  jet  center  line  downward  to  a 
point,  y  ,  where  the  mean  velocity  is  one-half  the  maximum  mean 
velocity  is  a  measure  of  the  jet  spreading.  The  lower  half  of  the 
jet  was  used  for  this  measurement,  in  order  to  reduce  the  influence 
of  the  free  surface  so  that  the  results  might  be  compared  with  the 
work  of  others.  Figure  7  shows  the  results  for  a  submergence  ratio 
of  H/D  =  6.9.  The  jet  expands  linearly  according  to 

y*/D  =  -  0.054  (x/D)  -  0.094  (2) 

(Correlation  Coefficient  =  0.994) 

This  equation  is  obtained  when  the  three  points  closest  to  the  jet 
source  are  included,  and  it  coincides  with  the  results  others  have 
reported  for  free  jets.  The  spreading  is  also  approximately  linear 
when  all  four  points  are  included  in  the  analysis.  The  equation  for 
this  is 
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y*/D  =  -0.134  (x/D)  +  0.728  (3) 

(Correlation  Coefficient  =  0.986) 

There  are  two  effects  that  make  this  situation  studied  different 
from  a  true  free  jet.  The  confinement  by  the  air-water  interface 
has  an  undetermined  effect  on  the  spreading  of  a  jet.  The  jet  was 
also  confined  from  below  by  the  solid  bottom  of  the  tank  which  was 
located  at  Hq/D  =  39.4  in  this  case.  The  point  farthest  downstream 
may  have  been  influenced  by  back  flow  along  the  tank  bottom,  though 
no  backflow  was  observed  directly.  Gutmark  and  Wygnanski( 1976 ) 
report  results  obtained  for  measurements  with  free  air  jets  in  air 
which  indicate  the  following  relationship. 

y*  -  0.1  (x/D  -  Xq/D)  (4) 

The  effect  of  the  free  surface  as  was  shown  above  is  to  distort  the 
centerline  of  the  jet  so  that  the  distance  from  the  centerline  to 
the  surface  becomes  smaller  as  x  increases.  This  could  account  for 
the  lower  coefficient  of  (x/D)  in  equation  (2). 

Turbulence  Quantities  Turbulent  velocity  fluctuations  were 
calculated  by  subtracting  the  mean  velocity  from  each  instantaneous 
velocity  value  in  the  time  series  which  were  calculated  from  the 
anemometer  output  data  as  described  above.  The  turbulent  velocity 
fluctuation  results  were  combined  in  the  appropriate  v;ays  and  time- 
averaged  to  give  the  Reynolds  stresses  and  turbulent  intensities. 

1)  Reynolds  stress 

The  Reynolds  stress  depends  on  axial  position  as  illustrated  in 
Figure  8  where  the  results  are  shown  for  a  submergence  ratio  of  6,9. 
Lines  drawn  in  the  figures  show  approximate  trends  and  are  not  "best 
fit"  lines.  These  dimensional  data  show  the  distortion  of  the  mean 
velocity  profile  due  to  the  presence  of  the  free  surface.  As  the 
jet  .^pxeads  out  downstream  from  the  source,  the  upper  and  lower 
parts  of  the  profile  become  less  and  less  antisymmetrical .  If  no 
free  surface  were  present,  one  expects  the  Reynolds  stress  profile 
to  be  antisymmetrical  about  the  centerline.  In  Figure  9,  the 
dimensionless  Reynolds  stress  is  plotted  as  a  function  of  the 
dimensionless  transverse  position.  This  shows  that  as  the  jet  moves 
downstream,  similarity  of  the  turbulence  is  approximately  maintained 
in  the  part  of  the  jet  below  the  centerline  away  from  the  free 
surface  but  not  maintained  in  the  part  of  the  jet  above  the 
centerline  toward  the  free  surface.  As  the  jet  moves  downstream, 
the  position  of  the  minimum  Reynolds  stress  moves  toward  the 
surface.  Based  on  the  Reynolds  stress  results,  one  cannot  say 
anything  definitive  about  the  position  of  the  centerline  of  the  jet. 
As  the  jet  moves  downstream,  the  magnitude  of  the  Reynolds  stress 
changes  first  increasing  and  then  decreasing. 

2)  Turbulence  Intensity 


The  turbulence  intensities  are  plotted  in  Figure  10 (axial 
turbulence  intensity)  and  Figure  11 (transverse  turbulence  intensity) 
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in  terms  of  dimensionless  position.  The  submergence  ratio  which 
applies  to  these  data  is  again  6.9.  As  the  jet  moves  downstream, 
the  dimensionless  intensity  profiles  are  reduced  in  size,  though  the 
magnitude  of  the  maxima  and  the  general  shape  are  approximately  the 
same  at  all  axial  positions  studied.  This  indicates  that  in  terms 
of  the  length  scale  y  similarity  of  the  turbulence  is  not  preserved 
in  the  planar  jet  flow.  Th.is  is  a  finding  different  from  that  based 
on  the  Reynolds  stress  results  for  which  similarity  was  observed  in 
the  lower  half  of  the  jet. 

The  uncertainty  is  rather  large  for  the  turbulence  intensities 
and  for  the  Reynolds  stresses  as  is  seen  from  the  data  scatter  of 
the  plotted  points.  The  tests  for  similarity  require  plotting 
dimensionless  values  which  are  calculated  by  multiplying  together 
two  velocity  fluctuation  values  which  in  turn  are  calculated  by 
subtracting  a  mean  velocity  from  a  measured  velocity — two  numbers 
which  differ  by  at  most  by  10  to  20  per  cent.  This  product  is  then 
dividt d  by  another  quantity  calculated  in  the  same  way.  The 

uncertainly  in  the  measured  velocities  is  approximately  3  to  15 

percent.  This  means  that  the  uncertainty  in  the  calculated 
turbulence  results  is  30  to  60  per  cent.  This  is  much  larger  than 

the  uncertainty  in  the  mean  velocities  because  it  is  likely  that  the 

uncertainties  associated  with  each  individual  velocity  measurement 
are  normally  distributed  and  that  the  time  averaging  process  keeps 
the  uncertainty  of  the  final  calculated  mean  velocity  approximately 
the  same  as  the  uncertainty  of  an  individual  measurement.  Since  a 
12-bit  A/D  converter  was  used  with  a  16-bit  microcomputer  for  data 
logging,  the  relative  uncertainty  associated  with  the  electronic 
data  processing  was  0.00024(1/4096).  All  the  uncertainty  in  the 
velocity  measurements  comes  from  sensor  calibration  and  from  not 
being  able  to  produce  a  true  stationary  jet  flow. 

3 )  Influence  of  the  Free  Surface 

From  the  intensity  results  reported,  it  is  not  possible  to  say 
anything  definite  about  the  effect  of  the  proximity  of  the  free 
surface  on  the  turbulence  intensity.  Others  have  postulated  that 
eddies  are  damped  by  the  free  surface  and  that,  compared  to  a  jet 
which  is  well-submerged,  the  transverse  intensity  should  be  reduced 
and  the  axial  intensity  should  be  increased.  Intensity  results  for 
three  different  submergence  ratios  are  shown  in  Figure  12.  The 
intensity  profiles  are  distorted  as  the  submergence  decreases.  The 
scatter  in  the  results  is  large  enough  to  mask  the  characteristic 
shape  with  two  maxima,  but  the  magnitude  of  the  intensity  values 
scaled  by  maximum  velocity  does  not  depend  on  submergence.  Figure 
13  shows  the  effect  of  submergence  on  the  Reynolds  stress.  The 
profiles  are  distorted  as  the  distance  to  the  free  surface 
decreases.  The  part  of  the  Reynolds  stress  profile  nearest  the  free 
surface  is  most  strongly  affected.  As  the  jet  is  submerged  farther 
from  the  free  surface,  the  characteristic  concave  outward  shape  of 
the  Reynolds  stress  profile  is  restored. 
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4)  Pressure  Fluctuations 

The  pressure  probe  proved  to  be  sensitive  enough  to  detect 
pressure  fluctuations  of  approximately  2000  Pascals.  This  is 
illustrated  in  Figure  14  where  the  difference  between  the  pressure 
at  the  sensor  tip  and  the  reference  pressure  is  plotted  as  a 
function  of  time.  Power  spectra  of  the  pressure  signals  are 
presented  in  Figure  15.  They  show  a  large  contribution  at  low 
frequencies  with  a  peak  at  32  Hertz.  When  the  probe  is  pointed 
directly  upstream,  this  peak  is  eliminated.  The  peak  is  greatly 
reduced  when  the  probe  is  pointed  directly  downstream.  We  could  not 
identify  any  noise  or  vibration  in  the  system  of  this  frequency,  but 
without  more  evidence,  we  should  not  infer  that  this  peak  indicates 
some  phenomenon  related  to  the  turbulence.  These  spectra  were 
originally  calculated  to  give  information  about  the  best  mounting 
angle  for  the  pressure  probe  to  avoid  an  effect  of  impact  in  the 
output  of  the  instrument.  We  found  that  the  best  position  for 
reducing  the  effect  of  impact  and  avoiding  the  effect  of  the  wake  of 
the  probe  support  was  135®  where  zero  degrees  indicates  that  the 
probe  is  pointed  directly  upstream. 

5 )  Autocorrelation  Functions 

Taylor's  hypothesis  allows  transformi’ig  spatial  and  temporal 
correlations  for  homogeneous  flows — those  which  have  constant  mean 
velocity.  The  autocorrelations — velocity /velocity  correlations  with 
one  of  the  velocities  displaced  in  time — were  calculated  using  the 
single-point  data  ,  so  that  we  could  determine  the  magnitude  of  the 
error  which  would  be  made  in  applying  Taylor’s  hypothesis  to  our 
non-homogeneous  jet  flow.  The  autocorrelations  show  the  expected 
shape  and  indicate  something  about  the  time  scales  associated  with 
the  turbulence.  The  details  of  the  autocorrelation  results  are 
discussed  in  the  attached  thesis  by  E.  Meaning — pages  35  to  41. 


C.  Two-Point  Results 

The  two-point  measurements  were  to  be  one  of  the  main  features 
of  the  whole  project.  The  information  is  essential  for  developing  a 
turbulence  closure  model  that  is  more  universal  than  single-point 
closure  models  such  as  the  k-epsilon  model.  Flow  features  in 
turbulence  which  persist  over  a  time  long  compared  to  the  time 
required  for  a  turbulent  fluctuation — coherent  structures — have  been 
identified  in  many  different  flow  fields.  Because  of  this,  it  makes 
sense  to  try  to  describe  turbulence  in  terms  of  the  range  of 
interactions  between  two  points  in  the  flow.  A  governing  equation 
for  the  transport  of  the  two-point  correlation  function  can  be 
written,  but  it  contains  triple  correlations  and  pressure-velocity 
correlations.  The  triple  correlations  must  be  modelled  for  all 
flows,  and  the  pressure-velocity  correlations  must  be  modelled  for 
non-isotropic  flows.  These  model  functions  must  be  evaluated  from 
experimental  data.  Once  the  model  functions  are  available,  they  are 
used  in  solving  the  differential  equation  for  the  two-point  velocity 
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correlation  which  in  turn  can  be  used  to  calculate  the  velocity 
field. 

One  objective  of  the  work  was  to  obtain  experimental  turbulent 
velocity  data  which  would  be  suitable  for  calculating  the  two-point 
pressure-velocity  correlations,  <puj^'>,  the  two-point  double 
velocity  correlation,  <u^u^'>,  and  the  two-point  triple  velocity 
correlations,  <Uj^u jUj^, '>.  These  measurements  were  not  successful, 
because  of  the  experxmental  uncertainties  which  affected  the 
results.  The  same  uncertainties  as  found  in  the  single-point 
measurements  were  present  with  the  two-point  measurements.  The 
uncertainty  introduced  by  subtracting  the  mean  value  from  all 
velocity  values  which  appeared  in  calculating  the  Reynolds  stress 
and  turbulent  intensity  also  appears  in  the  two-point  correlation 
calculations.  Using  two  hot-film  sensors  introduced  further 
uncertainty,  because  the  interactions  in  the  turbulence  took  place 
over  a  short  range  in  the  jet  flow.  The  wake  of  the  upstream  sensor 
seemed  to  interfere  with  the  measurements  obtained  from  the  down 
stream  sensor  when  the  two  sensors  were  placed  along  the  same  axis. 
These  difficulties  together  with  the  fact  that  the  relative 
uncertainty  is  increased  when  the  turbulent  velocities  are 
multiplied  together  in  calculating  the  correlations  combined  to 
decrease  the  signal-to-noise  ratio  to  the  point  where  no  definitive 
results  were  obtained  from  the  two-point  correlation  measurements. 

The  two-point  spatial  correlation  values  calculated  are  not 
accurate  enough  to  be  used  to  verify  the  results  of  closure 
modelling.  However,  the  results  are  consistent  with  measurements 
made  by  others.  The  complete  discussion  of  the  double  and  triple 
spatial  velocity  correlations  appears  on  pages  35  to  50  in  the 
attached  thesis. 


D.  Closure  Modelling  Work 

During  the  three  years  of  the  project,  the  theoretical  and 
modelling  work  has  been  focussed  on  developing  a  model  which 
describes  jet  turbulence  in  terms  of  parameters  which  depend  only  on 
Reynolds  number  and  are  therefore  fundamental  turbulence  parameters. 
To  develop  a  two-point  closure  for  anisotropic  turbulence,  it  is 
necessary  to  model  two  terms  in  the  two-point  velocity  correlation 
equation. 


<pu'>  related  to  anisotropy (  this  is  zero  for  iso¬ 
tropic  turbulence) 

<uuu'>  related  to  the  presence  of  scale  independent 
energy  transfer  modes  or  "preserved" 
structure 

These  quantities  may  be  modelled  in  terms  of  the  invariants  of  the 
separation  position  vector  and  a  unit  vector  in  the  direction  of  the 
flow.  Using  this  latter  vector  is  our  method  for  introducing 
anisotropy  into  the  model.  For  example,  the  pressure  velocity 
correlation  can  be  expressed  as — 
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<pujL’>  =  B3^(r  •  r,  —•  £*  k/  t)ki 

where  r  is  the  separation  vector  and  b  is  a  unit  vector  in  the 
direction  of  flow,  and 

P  =  P(X/  'ii'  “  '^i(2£  Et  t) 

As  an  approximation  for  the  model  being  developed,  let  there  be 
"small"  anisotropy.  This  allows  separating  the  model  into  two 
functions — an  isotropic  part  which  depends  on  the  separation  vector 
and  an  anisotropic  part  which  depends  on  the  preferred  direction, 

<puj^’>  =  Bj^{r  .  r,  t)rjL  +  B2(r»  b,  t)bj^ 

This  model  has  been  chosen,  because  it  provides  a  simple  way  to 
include  an  explicit  dependence  on  anisotropy  and  is  distinctly 
different  from  the  corresponding  expression  for  the  isotropic  case. 
It  has  been  determined  that  is  and  even  function  of  r  and  that  B2 
is  an  odd  function  of  b.  Because  of  mass  conservation,  Bj  and  B2 
are  not  independent  but  are  related  by 

Bi(r.r,  t)  =  -r-3|  r2  dr 

For  b  =  (1,0,0)  and  r  =  {r,0,0),  the  expression  for  the  pressure- 
velocity  correlation  gives — 

<p  u’>  =  B^  r  +  B2 

<p  v'>  =  <p  w'>  =  0 

The  triple  velocity  correlation  is  modelled  as  a  function  of  an 
energy  transfer  parameter,  k{r) .  It  has  been  determined  that  a 
reasonable  form  for  k  is — 

k{r)  =  a^  (r/v^  exp[  -a^Cr/t^  )2] 

where  is  a  "characteristic"  length  scale,  a2  defines  the  value 
of  {r/v^  )  for  which  energy  transfer  is  a  maximum,  and  a^  defines  the 
corresponding  amount  of  energy  transferred.  All  of  the  functions 
described  are  also  functions  of  time,  so  they  are  related  to  the 
decay  of  the  turbulence  a.^  the  water  moves  along  in  the  jet.  This 
decay  is  modelled  as  the  t>.  m  of  the  decay  of  an  isotropic  part  and 
the  decay  of  an  anisotropic  part  with  no  interaction  between  the  two 
parts . 

Because  of  the  experimental  difficulties  encountered  during  the 
course  of  the  project,  we  do  not  yet  have  any  suitable  data  for  the 
pressure-velocity  correlations  or  the  triple  velocity  correlations 
to  use  for  developing  the  closure  model  for  a  non-isotropic  flow. 

The  model  for  two-point  correlations  of  turbulent  velocity 
fluctuations  has  been  developed  using  the  two-point  triple  velocity 
correlations  for  isotropic  turbulence  reported  by  Stewart ( 1951 ) . 

This  theoretical  model  contains  two  parameters,  a  and  b,  which  are 
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pure  numbers.  Since  the  only  number  peculiar  to  the  analysis  of 
turbulence  is  the  Reynolds  number,  a  and  b  cannot  be  arbitrary 
parameters  but  must  depend  on  Reynolds  number  in  some  well- 
prescribed  way.  A  formal  representation  could  be  written 

a  ~  F{NRe)  and  b  -  G(NRe) 

F(...)  and  G{...)  may  be  different  for  different  turbulent  flow 
situations.  One  preliminary  result  of  fundamental  importance  is 
that  the  ratio,  a/b,  is  insensitive  to  changes  in  Reynolds  number. 
This  suggests  an  energy  cascade  autonomy  or  "self-governance"  which 
has  not  been  reported  in  the  literature.  The  properties  of  a  and  b 
are  as  follows. 

1 )  a  ^  0  and  b  0 

2)  a  <  b 

It  also  appears  that  {  b  -  a  )  is  a  measure  of  how  much  energy  is 
transferred  in  the  fluid  by  eddy  entrainment,  and  (b  +  a  )  is  a 
measure  of  how  much  energy  is  transferred  in  the  fluid  by  eddy 
deformation. 


E.  Publications  and  Conference  Presentations 


Besides  the  four  theses  listed  below  under  the  heading  of 
participants,  the  following  publications  have  appeared  as  a  result 
of  research  supported  under  this  project.  Copies  of  these 
publications  are  included  in  the  appendix. 

Trevino,  G.  1989  "Isotropic  Analysis  of  Grid  Turbulence", 

Int  Journal  Eng  Sci,  27,  pp.  1463-1471 

Trevino,  G,  1989  "On  the  Invariant  Functions  of  the  Turbulence 
Bispectrum",  Int  Journal  Eng  Sci,  27,  pp.  529-537 

Trevino,  G.  1989  "Energy  Transfer  in  Turbulence",  Physics  of 
Fluids  A,  1^,  pp.  2061-2064 

Trevino,  G.  In  Press  "On  the  Theory  of  Isotropic  Decay",  Appl  Mech 
Rev, 

The  following  papers  have  .~en  presented  at  technical  meetings 
and  conferences. 


Hubbard,  D.  W. (with  E.  J.  Hine  and  G.  Trevino),  "Experimental 

Techniques  for  Measuring  Turbulent  Velocity  Correlations  in  a 
Jet  F]ow  Near  an  Air-Water  Interface",  Paper  6AM6-3,  22nd 
Annual  Meeting,  Society  for  Engineering  Science,  Penn  State 
University (October  1985) 
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Hubbard,  D.  W.(with  E.  J.  Hine,  L.  R.  Heydenburg,  and  G.  Trevino), 
"Turbulence  in  Water  Jets:  A  Two-Point  Closure",  Paper  131e, 
1986  Annual  Meeting,  American  Institute  of  Chemical  Engineers, 
Miami  Beach (November  1986) 

Trevino,  G.  "The  Invariant  Functions  of  the  Turbulence  Bispectrum", 
International  Conference  on  Fluid  Mechanics,  Beijing, 
China(July  1987 

Hubbard,  D.  W.  (with  L.  R.  Heydenburg),  "Pressure-Velocity 

Correlation  Measurements  in  Turbulent  Liquid  Flows",  14th 
Annual  Mid-Western  Universities  Fluid  Mechanics  Retreat, 
Rochester,  Indiana (April  1987) 

Trevino,  G.  "On  the  Theory  of  Isotropic  Decay",  Pan  American 
Conference  of  Applied  Mechanics,  Rio  de  Janeiro, 

Brazil (January  1989) 
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G.  Plans  for  Future  Work 

Our  main  plan  for  future  work  for  this  jet  flow  turbulence 
project  is  to  improve  the  precision  of  the  experimental  data  so  that 
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the  closure  model  can  be  developed  for  non-isotropic  flows.  Wa  need 
to  control  the  flow  better  to  provide  truly  stationary  turbulence. 
The  jet  flow  we  used  for  our  measurements  was  produced  by  pumping 
water  into  a  reservoir  of  still  water.  Using  a  centrifugal  pump 
caused  small  fluctuations  in  the  flow  because  of  speed  variations 
due  to  minor  frequency  variations  in  the  line  power.  This 
difficulty  can  be  overcome  by  using  a  stand-pipe  water  supply  system 
which  is  available  in  the  Department  of  Civil  and  Environmental 
Engineering  at  Michigan  Tech.  This  stand-pipe  is  approximately  50 
feet  high  and  has  an  over  flow  weir  and  a  pumping  system  to  maintain 
a  constant  hydrostatic  head.  If  such  a  system  were  connected  to  the 
jet  flow  nozzle  unit,  the  upstream  pressure  would  be  constant  and 
there  should  be  no  fluctuations  in  the  mean  flow.  This  supply 
system  would  also  eliminate  any  difficulty  of  tem.perature  changes  in 
the  water  due  to  viscous  dissipation  when  the  nozzle  is  supplied 
directly  by  a  centrifugal  pump. 

The  turbulent  jet  is  supposed  to  be  two-dimensional,  but  there 
can  be  entrainment  at  the  sides  leading  to  transverse  velocity 
gradients.  A  refinement  which  would  eliminate  this  difficulty  would 
be  to  design  shields  for  the  sides  of  the  nozzle  so  that  water  would 
be  entrained  only  from  above  and  from  below. 

The  calibration  system  we  developed  during  the  course  of  this 
project  is  effective,  but  the  procedure  is  cumbersome  and  time 
consuming.  We  have  some  ideas  for  devising  an  in-situ  sensor 
calibration  device  which  would  be  a  moveable  laminar  jet.  Using 
this  kind  of  calibration  device  would  allow  leaving  the  sensor  in 
position  for  the  frequent  periodic  calibration  which  is  necessary 
when  making  measurements  in  water.  Calibration  data  would  be 
obtained  with  the  calibration  jet  in  one  position,  and  the  proper 
relationship  between  velocity  and  sensor  output  would  be  determined 
by  the  well-established  equations  giving  the  response  of  yawed 
sensor  elements. 

The  separation  distance  necessary  to  measure  non-zero  two-point 
velocity  correlation  functions  turned  out  to  be  quite  small.  This 
made  it  difficult  to  mount  the  sensors  directly  in  alignment  because 
of  the  arrangement  of  the  supporting  structures.  With  the  close 
spacing  needed,  the  wake  of  the  upstream  sensor  interferes  with  the 
signal  of  the  downstream  sensor.  This  makes  the  output  of  the 
downstream  sensor  noisy.  Both  the  alignment  and  proximity  problems 
could  be  eliminated  by  used  a  laser  doppler  anemometer (LDA)  as  the 
upstream  sensor.  This  is  a  nonintrusive  device  which  does  not 
disturb  the  flow.  The  downstream  sensor  could  be  a  crossed-sensor 
hot  film  probe  just  as  was  used  in  the  present  study.  This  idea  was 
first  suggested  by  Thomas  Sweam  of  the  Naval  Research  Laboratory. 
There  is  no  LDA  unit  available  at  Michigan  Tech. 

The  bleed  type  pressure  sensor  proved  to  be  sensitive  enough  to 
measure  the  pressure  fluctuations  which  occurred  in  the  jet  flow, 
but  the  device  was  not  rugged  enough  for  long  term  application  in 
liquid  flows.  If  the  pressure-velocity  correlation  measurements  are 
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continued/  a  more  rugged  pressure  sensor  is  needed.  A  piezoelectric 
pressure  sensor  might  be  appropriate. 

Computer  power  increases  all  the  time,  and  people  at  various 
laboratories  are  able  to  simulate  instantaneous  turbulent  velocity 
fields  by  solving  the  Kavier-Stokes  equations.  Results  are 
available  right  now  for  low  Reynolds  number  flows  in  simple  systems. 
It  would  be  interesting  to  try  to  extend  this  direct  simulation 
method  to  the  case  of  a  jet  flow  interacting  with  a  free  surface. 

If  that  effort  was  successful /  the  results  could  be  used  to  develop 
the  two-point  closure  scheme.  These  results  could  in  turn  be 
compared  with  experimental  data.  The  simulation  might  be  more 
effective  for  model  development /  because  the  velocity  values  are 
synoptic  and  closely  spaced  if  the  numerical  analysis  grid  is  fine 
enough . 
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Figure  1.  Jet  Geometry  and  Notation 


Velocity  (m/s)  Velocity  (m/s) 

Figure  2.  Effect  of  axial  position  on  mean 

velocity  profiles.  (D  =  6.35  mm,  H/D  -  6.9) 


Figure  3.  Dimensionless  mean  velocity  profiles 
(D  =  6.35,  H/D  =  6.9) 
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Figure  6.  Centerline  mean  velocities. 
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8.  Effect  of  axial  position  on  Reynolds 
stress  profiles.  (D  =  6.35  mm,  H/D  =-  6.9) 
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Figure  10.  Dimensionless  axial  turbulence  intensity 
profiles.  (D  =  6,35,  H/D  =  6.9) 
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Figure  12.  Effect  of  the  Free  Surface  on  Axial  and  Transverse 
Turbulent  Intensity  (D  »  12,7  ma,  x/D  »  11.0) 
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FIGURE15.  POWER  SPECTRUM  DENSITY 
AS  A  FUNCTION  OF  FREOUENCY 


APPENDICES 


A.  Copies  of  published  papers 

B.  M.S.  thesis  written  by  Eric  Menning 


Ml.  J.  Engng  Sci.  Vol  27.  No.  5.  pp  529-537.  1989 
Primed  in  Great  Bniain  All  rights  reserved 


0020-7225/89  $3  00  +  0  00 
Copyright  ©  1989  Pergamon  Press  pic 


ON  THE  INVARIANT  FUNCTIONS  OF  THE 
TURBULENCE  BISPECTRUM 

GEORGE  TREVINO 

Mechanical  Engineering-Engineering  Mechanics  Department.  Michigan  Technological  University. 

Houghton.  MI  49931.  U.S.A. 

Abstract — It  ts  established  that  the  tnvariant  functions  of  the  turbulence  bispectrum  can  themselves  be 
represented  in  terms  of  basis  functions  which  are  independent  of  time.  These  basis  functions  are  also 
scale-independent,  and  therefore  identical  for  large  eddies  and  small  eddies  altke.  Dynamical 
implications  of  this  result  are  discussed, 

'"What  we  observe  is  not  nature  itself  but  nature  exposed  to  our  method  of 
questioning"  [1]. 


1.  INTRODUCTION 

Several  papers  concerning  the  application  of  the  bispectrum  to  the  analysis  of  noise  and 
turbulence  have  been  presented  to  the  community  [2-13].  The  accepted  definition  of  the 
bispectrum  is  very  simple;  it  is  the  double  Fourier-transform  of  the  triple-moment  function,  and 
for  the  homogeneous  turbulence  is  typically  written  as 

5,;/(k„  k2,  /)  =  J  exp{-i(k,  •  q  +  k^  •  r)}5y,(q.  r,  t)  dq  dr,  (1) 

where 

5y/(q,  r,  t)  =  {u,(x,  r)u/x  +  q,  /)m,(x  +  r,  /)). 

The  bispectrum  is  the  three-point  generalization  of  the  power-spectrum  concept  already  fully 
developed  for  two-point  moments,  and  is  accordingly  an  energy-transfer  spectrum.  It  provides 
information  about  the  turbulent  flow  which  is  not  obtainable  through  the  customary 
energy-spectrum;  in  particular  it  supplies  knowledge  about  the  degree  of  nonlinearity.  The 
bispectral  tensor,  By;(ki,  k2,  r),  can  be  written  in  terms  of  the  invariant  functions  formed  from 
the  wave-number  vectors,  k,  and  k2,  and  for  isotropic  turbulence  characteristically  expresses  as 
(cf.  [14]) 

5y,(k„k2,f)=  S  A.,.F\-p(kM.  (2) 

Ai'*  1 

In  equation  (2)  the  f [,‘)''(k, ,  kj)  are  necessarily  real,  but  the  /l.v(- ■  •)  may  be  complex 
depending  on  the  nature  of  By;(ki,  k2, /),  cf.  eqn.  (1);  also  the  functional  form  of  each 
f‘,7'(ki,k2)  is  independent  of  the  corresponding  Asi-  •  ),  which  are  strictly  scalar  functions  of 
the  invariants,  kj  •  k,,  kj  •  k2,  k2  •  k2,  and  the  time  t.  In  fact  the  functional  form  of  each 
f{,7*(ki,k2)  is  determined  from  purely  invariant  theoretic  constraints  while  that  of  the  •  •) 
is  determined  strictly  from  fluid  dynamical  deliberations.  Proudman  and  Reid  [15]  have 
elegantly  shown  that  indice  symmetry  conditions  together  with  conservation  of  m?v<;  i educe  the 
complete  specification  of  5y,(k,,  k->,/)  to  two  (2)  independent  functions,  say  and 

>42 

The  purpose  of  this  note  is  to  establish  that  for  self-similar  decay  each  of  the  >4*^,  N  =  1,2, 
can  themselves  be  represented  in  terms  of  a  finite  number  of  basis  functions  which  arc  not  only 
scale-independent,  but  also  independent  of  time;  they  are  indeed  independent  of  initial 
conditions,  and  the  temporal  decay  does  not  at  all  affect  them.  This  result  suggests  that  during 
the  attendant  decay  there  are  some  aspects  of  energy-transfer  which  remain  fully  preserved 
(invariant),  and  are  therefore  in  some  sort  of  universal  equilibrium.  It  is  compatible  with  recent 
experimental  findings  that  “turbulence”  is  indeed  characterized  by  a  remarkable  degree  of 
structure  and  order  [16],  although  quantitative  evidence  of  the  herein  intended  structure 
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remains  yet  to  be  reported  in  the  literature.  This  “invariance  within  decay”  is  a  paradox  which 
rightfully  belongs  to  the  sub-structure  of  turbulence,  but  it  can  nonetheless  be  explained  in 
traditional  terms  provided  the  proper  formulism  is  posed.  Invariance  in  any  phenomenon  is 
important  because  it  fosters  the  search  for  those  properties  which  are  unique  to  the 
phenomenon  itself,  properties  which  are  independent  of  the  particular  choice  of  reference 
system;  it  asserts  that  the  substantive  laws  of  a  given  field  should  be  expressible  in  a  manner 
which  is  completely  independent  of  the  explicit  scales  chosen  to  measure  the  field  variables. 
For  turbulence,  which  characteristically  transfers  and  dissipates  energy  in  a  manner 
(seemingly)  consistent  with  the  scale-sizes  present  within  the  flow,  this  feature  is  fundamental, 
subtle,  and  new;  it  is  strikingly  “nonfluidic”  in  flavor.  From  a  dimensional  analysis  viewpoint, 
the  existence  of  such  structure  is  not  to  be  unexpected  since  energy-transfer  is  by  nature 
inherent  in  turbulence  decay,  but  length  scales  are  for  the  most  part  subjective  in  the  sense  that 
the  scientist  is  virtually  free  to  choose  for  analysis  any  particular  length  scale  of  practical 
interest.  Typical  length  scales  whose  growth  during  decay  is  commonly  studied  are  the  integral 
scale,  A.  defined  as 

A(/)  =  f  fir,  t)  dr, 

*'0 

fir,  t)  ~  longitudinal  correlation  function  of  turbulence,  and  the  dissipation  scale,  A,  defined  as 
A(r)=  although  other  equally  valid  length  scales  could  readily  be 

conceived  and  likewise  studied;  examples  of  such  are 

Z,,  =  A“'(r)[  r/(r,  r)dr 
■'0 

and  (  '1*^ 

L2  =  {A'‘(r)J^  rV(r,  r)drj  . 

Since  a  turbulence  Reynolds  number  can  be  formed  as 


Re 


/  Appropriate 
V  velocity  scale 


)( 


‘■'caf  )/ (''iscosily). 


the  proposed  result  infers  that  there  is  something  about  energy-transfer  in  a  self-similar 
turbulent  flow  which  is  independent  of  Reynolds  number,  regardless  of  how  the  flow  is  scaled. 
The  dimensional  analysis  approach  to  turbulence  invariably  assumes  that  a  turbulent  flow  can 
be  made  independent  of  Reynolds  number  only  if  it  is  “properly”  scaled.  For  the  case  of 
turbulence  generated  by,  say  a  wire  mesh,  the  present  implies  that  there  is  then  something 
about  energy-transfer  which  is  independent  of  grid  size  (gauge  invariance),  and  therefore 
identical  for  all  such  generated  flows.  From  an  analytical  dynamics  viewpoint  the  presence  of 
such  structure  is  also  not  to  be  unexpected  since  the  presence  of  sufficient  symmetry  in  a  given 
dynamical  system  almost  invariably  guarantees  related  constants  of  the  system  motion 
(conservation  laws).  The  knowledge  that  certain  functions  of  a  particular  system  actually 
remain  preserved  during  a  considered  motion  can  be  of  great  help  not  only  in  simplifying  the 
governing  equations,  but  also  in  eventually  leading  to  their  solution.  At  the  very  least  it  can 
lead  to  new  insights.  It  was  indeed  the  attendant  recognition  that  in  isotropic  decay  there  is 
inauspiciously  no  natural  length  scale  imposed  on  turbulence  which  prompted  Karinan  and 
Howarth  [17]  to  initiate  an  investigation  to  determine  if  closure  could  be  obtained  solely  from 
the  assumption  of  self-similarity  (self-preservance).  The  bispectrum,  therefore,  in  addition  to 
its  use  in  studying  such  mechanisms  as  spectral-transfer  and  vorticity  production,  appears  to 
also  be  a  tool  for  analytical  characterization  (and  study)  of  preserved  structure  during  the  decay 
of  self-similar  fluids.  Because  the  ensuing  has  its  genesis  in  topology-related  research  [18],  it  is 
accordingly  intended  only  as  a  description,  not  as  a  physical  theory,  of  turbulent  energy- 
transfer. 


2.  THE  CASE  FOR  SELF-SIMILARITY 


Since  the  assumption  of  self-similarity  is  vital  to  the  specific  results  of  this  communication,  it 
is  therefore  necessary  to  first  establish  under  what  conditions  such  behavior  would  be  expected 
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to  persist  in,  say,  a  grid-turbulence  or  even  a  jet-turbulence.  Historically,  the  assumption  of 
self-similarity  during  turbulence  decay  has  been  principally  a  mathematical  device,  employed 
solely  to  secure  definite  and  explicit  results.  Nonetheless  its  existence  and  subsequent  analysis 
can  be  equally  justified  from  purely  dynamical  system  grounds.  From  such  a  perspective, 
decaying  turbulence  is  a  time-varying  (nonlinear)  dynamical  system.  And  even  though 
time-varying  systems  can  be  regarded  as  “filters”  which  select  eigenfunctions  in  the  most  pure 
sense  rather  than  just  simple  sinusoids  (i.e.  they  discriminate  against  the  complex  exponentials 
as  the  natural  eigenfunctions  for  related  spectral  decompositions),  the  natural  tendency  m  their 
analysis  is  to  treat  them  as  systems  whose  properties  do  not  change  “appreciably”  with  respect 
to  time.  The  underlying  motivation  for  this  approach  is  that  classical  modal  analysis  is  not 
explicitly  restricted  to  systems  which  are  time-invariant  only;  it  is  only  necessary  that  the 
differential  equation  governing  the  system  as  well  as  the  boundary  conditions  be  separable  into 
a  function  of  time  alone  and  in  this  case  a  function  of  non-dimensionalized  length  alone,  i.e. 
that  temporal  changes  in  the  system  be  in  some  sense  “uniform”. 

For  purposes  of  mean-value  determination,  the  most  fundamental  mathematical  tool 
available  to  the  turbulence  practitioner  is  the  Reynolds  time-average.  In  any  real  flow 
inhomogeneity  is  always  present  and  therefore  the  best  estimate  to  the  mean  quantities  of  a 
steady-state  flow  are  provided  by  averages  such  as 

1 

Quu(Xi,X2)=  {u{Xut)l  u{Xi,  l)uiX2,  l)dt, 

1 

QuuuixuXz)^  t)u(x2,  0)  ir(x,.  f)u(jc:,  r)dr, 

etc.,  where,  in  practice,  measurements  of  flow  quantities  are  most  often  taken  along  the 
centerline  of  the  (spatial)  decay.  In  studies  of  turbulence  it  is  customary  to  introduce  a  velocity 
scale  defined  as  Ouix)  =  {Qu:,(x,  at)}*'*.  This  scale  permits  u{x,  t)  to  be  rewritten  as  u{x,  i)  = 
a{x}^{x,  t),  where  ^(x,  t)  is  a  normalized  turbulence  velocity  such  that  R{x,  x)  = 
(§(x,  t)^{x,  0)  =  1  for  all  X.  Note  that  o(x),  in  contradistinction  to  a  length  scale,  is  natural  to 
the  turbulence  itself  in  the  sense  that  a  particular  turbulence  will  automatically  define  its  own 
peculiar  <7(0:).  With  this  velocity  scale  in  hand,  Q„„(-  •  •)  and  0„„u(-  •  •)  immediately  rewrite  as 
X2)-*  C(x,  r)  =  a{x  -  rl2)o{x  +  rl2)R{x,  r)  and  Quuuixt.  Xj)— ►  C>(.r.  r)  =  (r(A:  - 
rl2)a{x  +  rl2)S{x,  r),  where  x  =  (x,  +X2)I2  and  r  =  X2-x^.  The  coordinates,  x  and  r.  can  be 
considered  as  uncoupled  principal  or  natural  coordinates  for  the  two  physical  coordinates,  Xi 
and  X2.  Unlike  modal  analysis  of  dynamic  systems,  however,-  the  principal  coordinates  here  are 
not  abstract  but  indeed  always  have  well-defined  physical  meaning.  Choosing  a  generic  length 
scale,  T](x),  it  is  further  customary  to  write  R(-  ■)-  R(r/rj),  meaning  that  the  algebraic  form 
of  R{-  ■  ■)  remains  unchanged  during  the  spatial  decay,  but  the  scaling  along  the  ab.scissa 
(r-a.xis)  has  stretched  according  to  the  growth  of  r]{-  ■  ■)  with  downstream  distance  from  the 
grid  (or  jet).  The  physical  condition  necessary  for  this  feature  to  manifest  is  that  energy  must 
cascade  in  such  a  manner  that  all  length  scales  of  the  turbulence  grow  with  x  at  approximately 
the  same  rate,  and  therefore  knowledge  of  the  growth  of  one  length  scale  e.ssentially  provides 
knowledge  of  the  growth  of  all  others:  accordingly,  there  is  no  accompanying  loss  of 
“information”  in  this  process  of  energy  transfer.  By  expanding  o(x  ±  r/2)  in  a  Taylor's  scries 
about  r  =  0,  the  result  is 

a{x  ±  r/2)  -  a(A-)  ±  ^  [  (r/2)  +  ■  ■  ■  ^  o(a){1  ±  rl2L  -!■■■■} 


L-/.(.v)  =  {a-(.v)k-}  }  ' 

I  I  or 

is  a  length  scale  over  which  a(-  •  ■)  may  be  assumed  to  vary  “appreciably".  Since  /?(■  ■  ■)  scales 
with  ?;(a),  the  net  effect  on  C(.r,  r)  of  the  different  length  scales,  L{x)  and  i]{x).  can  be 
characterized  by  writing  C(x.  r)  ~  o'ix){]  -  (r/2t])'{t]/l.)'}R{x,  r).  If  then  L»  r],  C(x,r)^ 
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a'{x)R{rlr])  for  those  values  of  r  for  which  r  -S-T).  Since 

lim  {R{x,  r))  =  0, 

r—*‘x 

if  L  is  prohibitively  larger  than  t],  C{x,  r)~  cr{x)R{rlr])  is  true  for  all  values  of  r.  Remote 
though  these  collective  conditions  are  from  anything  which  could  readily  be  established  in 
practice,  the  virtual  impossibility  of  their  manifestation  constitutes  no  objection  to  their 
empic>ment  in  a  purely  theoretical  investigation.  It  is  indeed  one  of  the  more  interesting 
features  of  any  physical  theory  that  it  invokes  idealized  behavior,  yet  the  eventual  results 
obtained  from  same  reveal  substantial  knowledge  about  real  phenomena.  Note  that  self¬ 
similarity  is  always  true  for  extremely  small  values  of  r.  The  formulation  is  ultimately  aligned 
with  isotropic  theory  by  expressing  downstream  distance,  ar,  as  x  -  Ui.  where  0  is  a  (constant) 
characteristic  mean-flow  velocity  (typically  a  centerline  velocity). 


3.  FORMULATION 

For  simplicity  consider  the  case  where  k2  =  ak,;  or -scalar.  For  this  case  the  A'"'’  become 
/l*''‘''(ki  •  k,,  k;  •  k2,  k,  •  k;,  r)~/4*‘'''’(A:,,  0-  A^  =  l,2,,  where  A:^  =  !k^l,  =  For  self- 

similar  decay  the  time-dependence  of  the  A"'’  is  carried  through  appropriate  scale-factors,  viz. 
o  ~  a(i)  and  rj  -  rj(r);  a  -  velocity  scale-factor  defined  through  a'(r)  =  J  i5y((ki,  kj,  r)  dk,  dk2, 
while  7j  -  subjectively  defined  length  scale  factor.  Accordingly,  A''^(A:,,/c2,r)~ 
r]k2)^  K2),  where  Kt  =  r]k,  and  K2  =  r]k2  are  time-dependent 

dimensionless  coordinates.  Proudman  and  Reid  (loc.  cit.)  have  established  that  ■  •)  must 
satisfy  certain  symmetries  in  the  arguments,  k^  and  k2.  viz.  that 

B,ji{ki,k2,t)  =  B,jiik2,k^,t)^  Bj,Xk2,  -k, -k2,  r).  This  symmetry  is  therefore  also  required  of 
the  scalar  functions,  k'2);  in  particular,  it  must  follow'  that 

K2)  =  f^i)  “  -  f^2)-  If  each  •  •)  is  expanded  in  a  Taylor’s 

series  about  *•,  =  1^2  =  0,  the  generic  result  is 

<p(fc,,  Jfj)  =  <P(0,  0)  +  0)'^!  +  0)/f;  +  ■■■  , 

where  ' ')  ' ')  denote  respective  appropriate  partial  derivatives.  Suppose  now 

that  the  Taylor's  approximation  is  restricted  to  the  zeroth-order  case,  viz.  the  case 
<p(/ri, /c-2)  ~  0)  “  constant.  It  is  obvious  that  <^>(0,0),  because  it  is  a  constant,  will 

automatically  satisfy  the  required  argument  symmetry.  Denote  (/)(0,  0)  by  IC2).  Suppose 

now  that  the  Taylor's  approximation  is  restricted  to  the  first-order  case,  viz  the  case 

K2)  =  K2)  +  0„.,(O,  0)*ri  +  ^,,,(0,  OjAT;- 

Since  //o(iCi,  k'2)  =  constant  already  satisfies  the  required  symmetry,  imposition  of  the 
Ki,  Ki-invertibility  on  <P(k,,  K2)  then  requires  that  (^)»-,(0,  0)  =  <^)k,(0,  0).  However,  imposition 
of  the  remaining  symmetry,  viz.  (p(K,,  K2)  =  0(k2, -k,  -  K2).  demands  that  <p„_(0.0)  =  0. 
Therefore,  there  is  no  linear  homogeneous  polynomial  in  Kj  ,  K2  that  under  any  circumstances 
satisfies  the  entire  required  k,,  Kj-symmetry.  Denote  {<^)»,(0,  0)k,  +  0»,(O,  0)K;}  by  H,(k,,  k-). 
Continuing  to  apply  the  same  procedure  to  the  second-order  approximation,  viz. 

(piKi,  K2)  =  W<,(k,,  K2)  +  /f,(K,,  k'2)  -(■  (l/2)(^,..,„(0-  0)^1 

+  ^/....-(O,  0)k,K2  -i-  (1/2)^, .,  .(0,  0)Ki 

produces  /^2(k,,  K2)  =  (any  constant)(K^  -t-  k,K2  +  K2).  The  procedure  applied  to  the  third-order 
approximation  yields  /^(k,  ,  k,)  =  (any  constant)(KtK;  +  KjKv).  Strikingly ,  though,  the  procedure 
applied  to  the  fourth-order  approximation  reveals  that  F/4(k,.  k'2)  =  (any  constant)(Kj -(- 
KiK2-t-K5)‘,  while  the  procedure  applied  to  the  fifth-order  approximation  supplies 
//5(k,,  K2)  =  (any  constant)(K7 -f  K1K2  +  kv)(k?K2 -t- k,K:),  i.e.  H4(- •  ■)  and  HX- ■  ■)  can  be 
constructed  from  H2i-  ■  •)  and  //,(•  •  •).  These  collective  results  are  not  as  surprising  as  they 
appear,  since  the  existence  of  a  finite  set  of  polynomial  basis  functions  for  any  sequence  of 
systematically  generated  algebraic  forms  was  established  in  an  elegant  theorem  by  David 
Hilbert  in  1890  [19j.  The  essence  of  the  theorem  formulated  by  him  is  that  if  an  infinite 
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sequence  of  algebraic  forms  (in  the  case  at  hand  the  sequence  is  //o(-  •  •),  //,(•  ■  •), 
Hzi'  ■•)>••)  constructed  according  to  some  general  over-riding  rule  (the  “rule”  here  is  the 
dominating  k^,  iCj-symmetry),  then  there  comes  a  point  in  the  construction  svhere  the  set  of 
available  independent  such  forms  is  exhausted  and  each  form  constructed  thereinafter  is  merely 
a  repetition/modification  of  what  has  already  preceded.  What  this  means  for  the  0‘"^(jc,,  /c,)  is 
that  K2)  must  be  constructable  from  some  combination  of  /c,)  and  K2), 

e.g.  •  •  •)  ~  Hli-  ••■)-!-  Hl(-  •  ■  ■),  while  //7(-  •  ■  •)  must  be  of  the  form,  W7(-  •  •  ■)  =  (some 
constant){//l(-  •  •  ■  ■  •)}•  No  other  fundamental  form  different  from  H2{-  ■  ■  •)  or  H,(-  •  •  •) 

may  occur  in  the  formulation. 


4.  DYNAMICAL  IMPLICATIONS 


D>ndmical  implications  of  the  foregoing  are  best  comprehended  by  examining  the  role  that 
the  bispectrum  pld>s  in  the  decay  of  total  turbulent  energy.  Keeping  in  mind  that  inertia  forces 
induce  energy  flow  between  different  wavenumbers  for  the  same  velocity  component,  it  is 
therefore  accordingly  necessary  to  analyze  only  the  simple  high  Reynolds  number  case 


du  1  du} 
3t  2  dx 


u~u{x,  t). 


The  two-point  correlation  method  duly  yields 


d{uu') 

dt 


+  -  («*u')}  =  0, 


where  u'  ~u{x  +  r,  t),  and  the  bispectrum  is  thus  introduced  through 


and 


ultimately. 


{uu'^)  =  j  B(ku  k2,  dk. 

{u^u')  =  j  B(ki,  k2,  f)e‘*-''d^,  dk^. 


S<P{k,  t) 

dt 


^  W{k,  t) 


where  t)  =  J  (uu'}e-‘*' dr  and  W(k,  t)  =  (ik/2)  f  (B(k.  k\ t)  -  B(k  -  k',  k',  t)}  dk'.  A 
dynamical  equation  for  the  bispectrum  (bispectral  decay  equation)  evolves  through 

du  d^u 

in  particular, 

°Bikuk2.t)  +  jt, A',  4.  kl)B(k, ,  it,.  0  =  0. 


An  immediate  feature  of  this  equation  is  the  explicit  presence  in  it  of  the  algebraic  form, 
//jfA'i,  It:)  -  (ik?  +  +  It:);  recall  that  the  necessary  presence  of  this  form  was  earlier 

established  on  purely  kinematical  grounds  yet  the  dynamics  of  the  phenomenon  establishes  its 
presence  independent  of  kinematics.  The  fact  that  //:(•  •  ■)  emerges  from  two  rather 
independent  lines  of  thinking  suggests  that  between  //:(■  •  ■)  and  //.,(•  •  ■),  //:(•  •  •)  plays  the 
stronger  role  in  the  transfer  of  energy  from  wave-number  k,  to  wave-number  k:  (or  vice  versa). 
The  solution  to  the  bispectral  decay  equation  is 

e(Ik,,  As,  /)  =  /3(Ik,,  As)e-='<*f**=*^-*=>' 

where  li(-  ■  ■)  is  determined  from  initial  conditions,  and  must  be  complex  (i.e.  have  both  real 
and  imaginary  parts).  This  solution  illustrates  an  interesting  but  not  unexpected  result,  viz.  that 
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(for  “low”  Reynolds  number)  the  decrease  of  energy-transfer  with  time  occurs  at  a  higher  rate 
for  the  eddies  with  the  larger  wave-numbers;  i.e.  energy-transfer  for  the  smallest  eddies  decays 
due  to  viscous  dissipation  much  faster  than  that  for  the  largest  ones.  Note  that  for  fc,  =  =  0, 

R(- ••)  s  constant.  Since  j3(- •  •)  must  also  be  expressible  in  terms  of  the  formulated 
fundamental  forms,  the  form  suggested  here  is 

.  .  fio  +  P,mk,  +  k,kiy 

■  /S;  +  -1- ^,^2  -1- ^2) 


This  form  substituted  into  the  bispectral  decay  equation  provides  the  right  A' -dependence  for 
•^{k,  /);  keep  in  mind  that  /3.v.  N  =  0,  I,  2,  may  all  be  complex. 

For  purposes  of  the  present  analysis,  however,  it  is  best  to  consider  B(-  ■  •)  as  a  function  of 
the  kinematic  invariant,  a  =  ^2/^1  =  K2/kx.  Accordingly, 


Biku  k2,  t)-* 


B{k,  a.  t)  = 


^0  +  Px(a  +  a^k” 

^2  +  +  O'  +  Cl!')k~  ^ 


2v(l  ♦o 


The  forms,  fii(cr)  =  1  -t-  ar  -1-  or*  and  h,i{a)  =  a  +  a',  can  be  thought  of  as  fully  invariant  (i.e. 
scale-independent,  time-independent,  and  Reynolds  number-independent)  “modes"  of  energy 
transfer;  note  that  there  are  only  two  such  modes.  In  each  mode  the  transfer  resembles  a 
simple  one  degree-of-freedom  system  (the  modes  being  discrete  and  independent  of  one 
another),  even  though  energy  is  transferred  through  a  continuum  of  wave-number  frequencies. 
Also  note  that  for  “small”  a,  h,{a) »  h„{a),  while  for  “large”  a,  h,{a)  =  hi,{a).  This  result  is 
consistent  with  the  widely  held  notion  that  the  effect  of  a  big  eddy  on  a  small  eddy  is  composed 
of  two  parts,  one  a  deformation  and  the  other  an  entrainment,  with  the  second  part  dominating 
when  the  eddies  in  question  differ  widely  in  size.  The  implications  here  are  that  h,(a)  - 
entrainment  while 

^  I  entrainment  for  “large”  cr, 

"  I  deformation  for  “small”  a. 


Furthermore,  for  “small  enough”  k,  i.e.  the  larger  scale-structure  of  the  turbulence. 


B(k,  or,  r)  = 


^2  +  Pih,(a)k‘ 


g-2%/n(o)C-f 


indicating  that  large-scale  energy-transfer  occurs  in  only  one  mode.  viz.  the  entrainment  mode. 
However,  since  entrainment  between  two  eddies  of  different  size  can  occur  only  when  the 
larger  eddy  totally  dominates  the  smaller,  not  only  do  the  “big  slow  eddies"  interact  very 
weakly  with  the  remainder  of  the  turbulence  (a  phenomenon  already  well-known  and 
understood),  but  the  present  establishes  that  they  also  interact  very  weakly  (if  at  all)  with  each 
other,  i.e.  eddies  of  roughly  the  same  size  cannot  entrain  one  another.  This  is  perhaps  what 
accounts  for  the  fact  that  the  turbulence  spectrum,  ^{k,  t),  has  a  somewhat  permanent  form 
around  k-0.  Because  the  mode,  h„{a),  does  not  participate  in  the  large-scale  transfer  of 
energy,  the  dynamical  system  (i.e.  the  turbulent  fluid),  for  this  wavenumber  domain  behaves 
like  a  degenerate  or  semi-definite  system;  from  this  perspective  h„{a)  is  like  an  ignorable 
coordinate  and  the  corresponding  “constant  of  the  motion”  is  the  permanence  of  the 
large-scale  turbulence  structure.  For  the  larger  values  of  k,  however. 


Bik,  a,  0  = 

Pikiia) 

indicating  that  for  the  smaller  scale-structure  both  modes,  one  mode  (viz.  h„{a))  determined 
from  pure  kinematics  only,  participate  in  the  transfer.  This  condition  persists  even  for  the 
smallest  scales.  The  different  modal  behavior  for  small  and  large  scales  reveals  that  even 
though  kinematic  similarity  is  assumed,  formulated,  and  pursued,  the  system  nonetheless 
refuses  to  behave  dynamically  similar.  The  presence  of  the  mode.  in  the  respective 

expressions  for  large-  and  small-scale  energy-transfer  alike  supports  the  assertion  made  earlier 
that  this  mode  plays  the  stronger  role  in  the  transfer  of  energy. 
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5.  DEEPER  IMPLICATIONS 

Since  self-similar  turbulence  must  decay  in  a  fashion  which  preserves  the  scale-  and 
time-independent  forms,  /i/(a)  and  hii{a)  =  h,(a)  -  1,  and  since  the  bispectrum  for  any  such 
decay  must  be  algebraically  expressible  in  terms  of  h,(a)  and  h„{a),  then  all  universal 
ramifications  of  same  are  embodied  in  the  single  absolutely  invariant  function,  h,{a),  and 
eventually  in  the  invariant,  a.  The  physical  implication  of  this  result  is  best  comprehended  by 
recalling  a  simple  experiment  in  physics  which  culminates  in  a  precise  operational  definition  of 
the  mass  of  a  body.  The  experiment  involves  placing  a  spring  between  two  isolated  masses,  A 
and  B,  pushing  them  together  so  that  the  spring  is  compressed  (the  spring  is  not  rigidly 
attached  to  either  of  the  masses),  and  then  releasing  them.  By  repeating  the  experiment  a 
number  of  times,  each  time  with  a  spring  of  different  constant,  it  will  be  found  that  although 
the  recoil  velocities  of  A  and  B  will  be  different  each  time,  the  ratio,  v^/vg,  will  be  identical 
for  all  cases,  i.e.  v^/vg  =  invariant,  regardless  of  the  amount  of  initial  energy  in  the  system. 
The  fact  that  this  result  provides  a  means  for  unambiguously  defining  the  mass  of  a  body,  viz. 

-  vjvg,  and  that  it  is  consistent  with  the  principle  of  linear  momentum  conservation, 
is  well  established.  Put  simply,  the  invariance  of  v^lvg  means  that  energy  cannot  be  transferred 
from  the  potential  energy  of  the  spring  to  the  kinetic  energy  of  the  masses  in  an  arbitrary 
fashion  but  must  indeed  be  transferred  in  a  manner  consistent  with  linear  momentum 
conservation.  In  thermodynamics,  a  correspondingly  simple  experiment  will  establish  that  in  a 
Carnot  cycle  the  heat,  Q^,  flowing  in  at  a  high  temperature,  T^,  is  related  to  the  heat,  Q^, 
flowing  out  at  a  lower  temperature,  7),,  by  the  relation 

~  =  ^  ”  Invariant: 

Ql  Tl 

the  invariance  of  this  ratio  is  a  consequence  of  the  principle  of  conservation  of  entropy,  i.e, 
heat  is  not  transferred  arabitrarily  through  the  engine.  In  geometrical  optics  the  fact  that  the 
sines  of  the  angles  of  incidence  and  refraction  stand  in  a  constant  ratio  to  each  other,  regardless 
of  the  magnitude  of  each,  leads  to  the  formulation  of  Snell's  Law.  In  the  study  of  relativity, 
Mermin  [20]  has  invoked  the  invariance  of  the  ratio  of  collinear  lengths  for  both  the  fixed  and 
moving  observers  to  cleverly  derive  the  Lorentz  transformation  equations  without  making  use 
of  the  principle  of  the  constancy  of  the  velocity  of  light.  It  appears,  therefore,  that  the 
invariance  of  a  =  k2lku  as  herein  formulated,  is  a  manifestation  of  some  conservation  law  or 
principle  which  governs  the  transfer  of  energy  in  a  decaying  turbulence.  Such  a  conservation 
law  must  describe,  for  self-similarity,  the  means  of  (scale-independent)  “communication” 
between  turbulence  events  of  different  scale-sizes  such  that  in  the  limit,  events  distinguished  by 
a  great  disparity  in  scale-size  should  exert  minimal  influence  on  one  another.  Moreover,  such  a 
principle  is  consistent  with  the  traditional  wisdom  that  turbulence  does  not  transfer  energy 
arbitrarily  but  indeed  does  so  in  some  well-prescribed  fashion.  If  we  think  of  a  “parcel”  of 
turbulence  of  a  particular  scale-size  as  a  dynamical  system  which  receives  energy  from  those 
parcels  of  turbulence  of  larger  scale-sizes,  dissipates  some  of  this  energy,  and  transfers  the  rest 
to  those  parcels  of  smaller  scale-size,  then  a  possible  such  law  is 

E^R{k\,  ^2)  ^kj  \kJ 

in  this  expression  £r(*i)  is  the  total  energy  transferred  by  a  turbulence  parcel  of  scale-size, 
kr‘.  to  all  other  turbulence  parcels  of  smaller  scale-size,  E^ik^  k,)  is  that  portion  of  £r(^i) 
which  is  eventually  received  by  a  pracel  of  scale-size,  kj y  is  a  non-adjustable  constant  greater 
than  zero,  and  k2>ki  (see  Fig.  1).  The  /-dependence  of  both  £r(’ '  )  and  £r('  •  •)  is  explicitly 
deleted  in  the  above  expression  in  order  to  emphasize  the  time-independence  of  the  ratio, 
Et{-  ■  ■)/Eg{-  ■  ■).  The  ultimate  consequences  of  such  a  law  (or  others  like  it)  need  to  be  fully 
investigated. 
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Fig.  1.  D>namical  system  schematic  of  a  turbulence  parcel  of  scale-size.  ki'.  dissipating  energy  and 
transferring  energy  to  turbulence  parcels  of  smaller  scale-sizes. 


6.  CONCLUDING  REMARKS 

In  closing  it  remains  to  affirm  that  regardless  of  the  nature  of  the  decay — whether  self-similar 
or  non-self-similar,  isotropic  or  non-isotropic — some  aspect(s)  of  the  turbulence  will  always  be 
preserved.  This  follows  from  the  nature  of  scaling  effects  (transformations)  in  general  that 
regardless  of  how  a  geometric  figure  may  be  deformed,  something  about  it  nevertheless 
remains  intact  and  unchanged  during  the  deformation  process  [21].  The  geometric  figure  in 
question  here  is  the  shape  of  the  longitudinal  correlation  function  f{r,  t),  of  turbulence,  and 
how  this  function  evolves  in  time  determines  where  the  first  evidence  of  preserved  structure  can 
be  found.  For  self-similar  decay,  viz.  f{r,  t)  ~f{rlr)),  where  t)  ~  r?(/)  only,  invariant  structure, 
as  herein  formulated,  is  found  in  the  energy-transfer  mechanism.  This  invariance  accordingly 
admits  a  representation  of  this  mechanism  in  terms  of  a  finite  number  of  transfer  modes  which 
are  independent  of  the  magnitudes  of  the  coordinates  originally  chosen  to  describe  the  motion. 
In  each  mode  the  transfer  resembles  a  one-degree-of-freedom  system  (with  no  interplay 
between  the  modes),  even  though  energy  is  transferred  and  dissipated  in  the  entire  system 
through  a  continuum  of  wavenumber  frequencies.  For  a  more  general  strain  of  decay,  say 
where  possibly  tj  ~  7j(r,  r),  preserved  structure  will  necessarily  be  found  only  in  the 
higher-order  moments  [18].  This  means  that  independent  of  adventitious  disturbances, 
preserved  structure  of  some  sort  will  exist  in  any  and  all  types  of  decay;  furthermore,  it  can  be 
definitely  concluded  that  the  degree  of  particularly  endowed  structure  is  not  affected  by  the 
“age”  of  the  flow. 
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Abstract — The  principles  of  self-similar  isotropic  turbulence  are  applied  to  the  analysis  of  grid-decay 
In  particular  it  is  theoretically  established  that  energy-transfer  occurs  in  only  tsso  distinct  invanant 
modes.  The  decay  constant  of  third-order  correlations  is  posed  in  terms  of  the  '■amounts”  of  each 
mode  present  in  the  decay.  This  constant  is  determined  empirically,  and  is  found  to  be  relatisely 
insensitive  to  changes  in  Re  The  related  decay  of  turbulence  intensity  is  computed. 


1.  INTRODUCTION 

The  theory  of  isotropic  turbulence  centers  upon  solving  the  dynamical  equation. 


\or  r'  \ar‘  ror! 


to  explicitly  obtain  the  time  dependence  of  a  peculiar  length  scale.  7/(f),  in  terms  of  the 
time-decay  of  the  velocity  scale.  a(t).  The  velocity  scale  is  a  parameter  which  is  “natural"  to 
the  turbulence  itself  (typically  a  turbulence  rms  value)  while  r/(/),  on  the  other  hand,  is  a 
“subjective"  parameter  in  the  sense  that  the  analyst  is  virtually  free  to  choose  for  analysis  any 
length  scale  of  particular  interest.  Typical  length  scales  whose  growth  during  decay  is  studied  in 
the  literature  are  the  dissipation  scale,  defined  as 


or  the  integral  scale,  defined  as 


5-/1"'- 

or 

Ml)  =  I  f{r.  t)  dr. 

At 


The  only  natural  length  scale  of  turbulence  is  characteristically  that  imposed  by  the  mechanism 
which  generates  the  turbulence,  say  for  example,  the  mesh  size  in  a  grid-turbulence  or  the  jet 
diameter  in  a  jet-turbulence,  but  this  scale  is  (almost)  always  constant  and  essentially 
inappropriate  for  decay  analyses.  In  a  practical  scenario,  o(- •  ■)  is  invariably  approximated 
from  steadv-state  turbulence  measurements  as 


(x.  I)  dr. 


where  x  is  the  downstream  distance  from  the  turbulence-generating  mechanism:  further, 
X  -  Ot.  0  =  convectivc  mean-flow  velocity.  Similarly,/!-  •  ■)  and  k{-  ■  ■)  are  approximated  as 


1  r' 

o-{x)f{r.  .v)  -=  -  J  M(.tc.  t)u{x  -s  r.  r)  dr 
1 

a\x)k(r,  .rr=  -  I  i)Mx  +  r.  t)  dr: 


in  all  these  approximations,  u  is  the  .v -component  of  the  three-dimensional  turbulence  velocity, 
u  =  («.  i'.  »v).  T  is  a  long  averaging  time,  and  all  measurements  arc  typical'y  taken  along  the 
flow  centerline  where  because  of  the  attendant  zero  mean-flow  gradients  the  isotropic  theory  is 
most  applicable. 
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In  order  to  solve  equation  (1)  from  a  purelv  theoretical  perspective,  some  rather  ad  hoc  but 
nonetheless  plausible  assumption  needs  to  be  made  about  the  function,  k  k[r.  x)  -k(r.i). 
.An  especially  useful  one  which  has  on  occasion  been  invoked  (Ij  is  that  of  self  similantv  in 
both  /(r.  /)  and  kir.  r).  historically,  this  assumption  requires  that  /(•  -)  and  k(-  -)  b<iih  scale 
in  the  same  way  with  a  length  scale  which  is  a  function  of  time  alone  (and  completely 
independent  of  r).  but  it  will  be  shown  momentarily  that  such  a  requirement  is  sirictlv 
unnecessary  and  /(-  •  -)  and  k{-  ■  )  may  in  fact  each  scale  with  totally  different  length  scales 
without  this  condition  affecting  the  fluid  decay.  Fividence  that  /(•■•)  and  k{- ■  in 

grid-turbulence,  do  indeed  scale  in  different  ways  has  already  been  provided  by  Stewart  (2) 
when  he  plotted  values  of  h{-  -  -)  vs.  r/.V/.  M  ~  mesh  sue;  Stewart's  graphic  description  of  the 
large-scale  behavior  of  h{r/M)  is  that  it  ".  .  .  falls  remarkably  slowly  toward  zero,  and  does  not 
disappear  until  well  beyond  the  position  at  which  the  double  correlation  becomes  undctec- 
ta'ole."  Since  A{-  ■  ■)=  -2h(-  ■  ■)  the  same  description  is  therefore  applicable  to  k{-  -  -).  The 
decay  equation  which  emerges  from  the  self-similarity  assumption  is 


d(a-C,r/)  ,  , 


(2) 


where  Q  is  a  constant  which  is  unity  if  t]  =  and  greater  than  unity  if  rj  A.  C,  is  a  constant 
defined  as 

and 

Ct=4[  r*’A{r, /) dr. 

Jli 

Note  that  for  self-similarity,  k(r.t)-  k(rtL).  L-any  time-dependent  only  length  scale,  and 
accordingly 

f  r-‘A(r'Z.)dr=  f  [rlL)‘'k{riL)d{rlL).  (3) 

.Also  note  that 


is  not  independent  of  C,.  The  subsequent  decay  law  is  found  by  assuming  rj(r)-f  and 
a(r)-r".  from  which  it  follows  that  for  self-similarity  the  only  possible  values  of  r  and  d  are 
c=0.5=  -d.  This  same  result  was  originally  established  by  Dryden  j.'].  .-X  peculiar  conse¬ 
quence  of  this  result  is  that  the  turbulence  Re.  defined  as  Re  -  at]  v.  is  constant  throughout 
the  decay. 

On  occasion,  equation  (1)  is  addressed  in  its  inviscid  (infinite  Re)  form.  vu. 

Bsll 

3i 

yielding  for  the  self-similar  c.oe  the  family  of  solutions,  c  ~  \  ~  d.  Re  -0(r' '•■').  However, 
such  a  family  of  solutions  should  be  viewed  with  cauiuvn.  since  it  is  mctmsistent  wish  the 
fund.imcntal  precept  of  turbulence  that  mass  conservation  tv'gcthcr  with  momentum  cv-nserva- 
tion  duly  provides  four  cquatum.s  and  four  unknowns,  thus  rigorously  requiring  one  solution 
only  to  the  decay  problem.  Furthermore,  since  for  "small"  r.  k{  ■  )  -  O(r').  any  Mvlution  to 
the  inviscid  equation  cannot  be  valid  for  the  small-scale  structure  of  the  turbulence  (and 
therefore  not  a  viable  solution  at  all);  recall  that  for  small  r.  f{r.  r)  1  -  (r  rj)', 

Becau!.c  the  theory  of  turbulence  strictly  prohibits  the  presence  in  the  analysis  of  any 
adjustable  or  "tuneable"  con'tants.  a  complete  solution  of  the  dynamical  equation  is  therefore 
not  achieved  until  a  value  for  is  determined,  or  at  least  a  viable  method  for  so  determining  it 
explicitly  formulated.  In  other  words,  since  the  only  flow-rel.ited  number  unique  to  turbulence 
is  Re,  any  other  flow -related  numbcr(s)  present  in  a  given  turbulence  analysis  necessarily 
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has(have)  to  ultimately  depend  through  some  well-prescribed  method  on  Re.  ITie  fomulation 
of  such  a  method  and  its  implications  for  grid-turbulcncc  are  presented  in  the  following.  The 
fact  that  Q  is  formed  from  an  integral  of  k(-  ■  •)  fortunately  suggests  that  aside  from  some 
fundamental  limiting  constraints,  viz. 

lim  {/c(-  ■  )}  -Ofr'f 
and 

lim  {k{-  •  •)}  =  0. 

r— ** 

the  exact  functional  form  of  k{-  ■  •),  although  important  is  not  devastating.  Since. 

Q-i  r-A:(---)dr. 

•'0 

k(-  ■  ■)  also  does  not  have  to  be  integrable  over  all  positive  values  of  r;  i,e.  because  k{-  ■  •)  is  an 
odd  function. 

I  /c(---)dr  =  0 

J  -X 

while 

f  k(-  -)  dr 

■'0 

itself  may  be  undefined.  Furthermore,  because  of  the  rote  it  plays  in  decay, 

f  r-'k{‘--)dr 

must  always  be  less  than  zero. 


2.  FORMULATION 


The  transfer  of  energy  in  a  system  obeying  nonlinear  equations  is  a  dynamical  process 
explicitly  dependent  upon  the  equations.  Bounds  on  the  transfer,  however,  can  be  imposed 
from  pure  kinematics  only,  although  not  much  more  can  be  said  about  the  functional  form  of 
the  transfer  without  employing  dynamical  arguments.  The  first  step  in  the  formulation  is 
therefore  to  invoke  kinematics  to  establish  said  bounds  and  accordingly  develop  a  repre¬ 
sentative  algebraic  form  for  k{r,  i)~k(r/T}):  recall  that  the  exact  algebraic  form  of  k{-  ■  ■)  is 
not  central  to  the  problem.  To  fundamentally  accomplish  this  while  simultaneously  saying  as 
little  as  possible  about  the  transfer  mechanism  itself,  it  is  here  proposed  to  kinematically 
describe  the  three-point  correlation. 


Q{q-  r,x) 


)i/(x  -r  q.  t)u{x  +  r,  t)  di. 


where  again  i)  is  only  one  component  of  the  steady-state  three-dimensional  turbulence 
velocity  field  and  x  a  one-dimensional  position  vector.  In  a  strict  nonhomogeneous  turbulence, 
such  as  that  generated  by  (say)  a  grid,  the  correlation.  Q{q,  r,  x),  accordingly  scales  as 

Q{q.  r,  x)  =  o(x)0{x  +  q)o{x  -i-  r)s{qlq,  r/q). 

However,  the  self-similar  theory  assumes  that  the  length  scale,  /.  which  scales  a(-  •  ■)  is  “large" 
in  comparison  to  ij;  large  enough  that 


< 

a(x)o(x  -i-  q)o{x  +  r)s{qh],  rjq)  -  (7(x)a|x 


xsiqh],  rlq)=-  o\x)s(qh},  rfq). 

Through  its  dependvtnce  on  q  and  Q(-  •  -)  analytically  described  how  a  turbulence  eddy  of 
scale-size,  q,  "inteiacts  with”  (transfers  energy  to)  a  turbulence  eddy  of  scale-size,  r,  or 
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vice-versa;  i.e.  Q(-  ■  •)  quantities  how  turbulence  eddies  of  differing  size  “communicate”  with 
one  another  and  thus  measures  the  “strength”  of  the  nonlinearity.  For  this  reason,  Q{-  ■  ■) 
must  approach  zero  as  the  difference  between  q  and  r  grows  without  bound. 

A  simple  property  of  Q{q,r,x)  is  (cf.  Ref.  [4])  Qiq.  r,  x)  =  Q{r,  q,  x)  =  Q{-q,  r  -  q,  x)\ 
this  r-argument  symmetry  must  therefore  also  be  satisfied  by  s(- ■  ■).  If  5(- ■  •)  is  then 
expanded  it.  '  Taylor’s  series  in  q  and  r,  the  formal  result  is 

s{q/q,  rlq)=s{0,  0)  +a,g  +  b^r  +  +  C;^r  +  ■  •  • . 

where  the  a„,  b„,  c„  n  =  1,2,3,  ,  represent  the  appropriately  evaluated  partial  derivatives 

of  5(-  •  •)■  Following  the  procedure  outlined  in  Refs  [5]  and  [6],  it  can  be  readily  shown  that  due 
to  the  required  argument  symmetry  in  s(-  •  •),  aside  from  ^(0,  0)  which  automatically  satisfies 
the  required  symmetry,  the  only  algebraic  forms  which  can  appear  in  the  expansion  are  the 
homogeneous  polynomials, 

F\{q,  r)  =  q- -  qr  +  r  and  Fi{q,  r)=2q^ -3q~r -3qr- +  lr^, 

and  those  higher-order  homogeneous  polynomials  which  can  be  formed  from  f,(-  •  ■)  and 
^2(’ •  •).  c-g-  f^i(' ■  •).  ' ')h{' ■ ')<  Fl{- ■  ■),  etc.  No  other  fundamental  algebraic  form 

different  from  fi(-  •  •)  and  F2{-  •  •)  satisfies  the  required  symmetry  and  can  therefore  not  appear 
in  the  expansion  of  s(- ■  ■)■,  specifically,  the  (^,  r)-dependence  of  s{qlq,_rlq)  can  be  fully 
characterized  by  F|(- ■  •)  and  F2(- •  •),  i.e.  as  far  as  q  and  r  only  are  concerned, 

5{F,(9,  r),  Fjiq,  r)}. 

The  physical  significance  of  F,(-  •  •)  and  F2(-  ■  •)  is  determined  by  writing  Fi{q,  r)  =  (a"  -  a  + 
l)r^  =  hi{a)r  and  F2(^,  r)  =  2(o''^  -  l.Sa"  -  1.5q  +  l)r^  =  2/i2(ar)r\  where  a  =  qlr  is  a  scale- 
independent  parameter  which  is  an  invariant  of  the  decay,  i.e.  during  the  decay  a  =  qlr~* 
i,qlq)l{rlq)  =  a,  or  in  words,  as  the  big  eddies  are  “scaled  down”  into  little  eddies  in  a 
unilaterally  self-similar  fashion,  the  ratio  of  two  eddies  remains  preserved  and  unaffected. 
Consequently,  the  forms  hi(a)  and  hzia)  are  also  invariants  of  the  decay,  and  upon  further 
reflection  are  seen  to  be  identical  for  both  homogeneous  turbulence  and  nonhomogeneous 
turbulence  alike.  It  is  instructive  to  momentarily  divert  here  to  point  out  that  in  a  1981 
communication.  Deissler  [7]  concluded  that  even  for  a  general  nonhomogeneous  turbulence, 
certain  terms  in  the  two-point  stochastic-averaged  Navier-Stokes  equation  can  be  interpreted 
as  “transfer”  terms.  He  accomplished  this  by  introducing  a  third  point,  defined  as  x„  = 
X  F  n{x'  -at),  Os/1  <  1,  into  the  two-point  geometry  of  the  flow  and  cleverly  deducing  that  the 
appropriately  formulated  “turbulence  self-interaction”  terms  produce  zero  total  contribution  to 
the  decay  of  (mm');  u~u{x,  t),  u'  ~  u{x',  t),  x'  =  x  +  r.  The  impact  of  his  result  (duly  pointed 
out  by  him)  is  that  turbulence  self-interaction  plays  the  same  role  in  nonhomogeneous 
turbulence  that  it  does  in  homogeneous  turbulence,  viz.  to  transfer  energy  from  one  part  of 
wave-number  space  to  another  without  changing  the  total  amount.  Furthermore  he  makes 
perfectly  clear  that  the  interpretation  of  self-interaction  as  energy-transfer  is  correspondingly 
conceptually  independent  of  the  condition  of  homogeneity;  recall  that  the  condition  of 
homogeneity  is  typically  invoked  in  making  the  “transfer”  interpretation.  More  fundamentally, 
though,  Deissler's  result  suggests  that  the  two  types  of  turbulence,  through  the  energy-transfer 
mechanism,  are  somehow  “unified”.  In  the  present,  the  forms  h,(a)  and  h^ia),  because  of 
their  insensitivity  to  the  condition  of  homogeneity/nonhomogeneity,  are  collectively  one 
manifestation  of  this  unity.  Because  the  entire  transfer  mechanism  can  be  characterized  in 
terms  of  these  time-independent,  scale-independent  forms,  these  forms  then  can  be  thought  of 
as  fully  invariant  dynamical  “modes”  of  energy-transfer.  Note  that  there  are  only  two  such 
modes,  and  that  the  “two-mode”  behavior  is  a  direct  consequence  of  argument  symmetry  only 
in  ^(-  •  •),  while  the  formulated  scale-independence  is  a  consequence  of  the  imposed 
self-similarity.  Accordingly,  it  is  not  at  all  unreasonable  to  envision  the  herein  considered 
transfer-modes  analogous  to  the  modes  of  classical  heat-transfer,  viz.  conduction,  convection, 
and  radiation,  save  that  here  there  are  two  distinct  transfer  modes  only  in  contradistinctiuii  io 
the  three  of  heat-transfer.  The  two-mode  behavior  is  consistent  with  the  already  known  [8] 
phenomenon  that  the  effect  of  a  large  eddy  on  a  smaller  eddy  is  composed  of  two  parts — one  a 
deformation,  the  other  an  entrainment.  As  in  classical  modal  analysis.  /i|(a')  anJ  /i2(q')  are 
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Fig.  1.  Schematic  of  invariant  energy-transfer  modes 


intended  to  remain  intact  and  unaltered  during  the  decay  and  an  orthogonality  condition 
between  them  is  found  by  setting  p  ~  a-  0.5,  whereupon  ♦  /i(j3)  =  /3"  +  0.75, 

m  -  2.25)3,  and 

f"  /.(/3)/>(/3)d/3  =  0  foranye. 

J-e 

Plots  of  /tiCo")  and  /12(a)  for  0^  a{=qlr):^  1  are  shown  in  Fig.  1.  From  purely  dimensional 
arguments,  then  it  follows  that  s{-  ■  •)  =s{3{t)hi{a)r^,2Ait)h2ia)r^)  where  A(t)  =  a^lr]^, 
B{t)  =  b~!r\^,  {a,  b)~flow  constants  which  depend  only  on  Re,  and  thus  on  initial  conditions, 
and  the  e.xplicit  algebraic  expression  for  s{-  •  •}  varies  from  flow  to  flow.  i.e.  5{-  •  will  have  a 
given  form  for  grid-turbulence,  another  form  for  jet-turbulence,  etc. 

It  is  a  rather  simple  analysis  to  deduce  that  the  flow  constants,  a  and  b,  define  explicitly  “how 
much"  of  each  transfer-mode  is  present  in  the  decay.  Since  for  self-similarity  Re  =  constant,  it 
follows  that  a  =  constant  =  fo(R6)>  ^  =  constant  =  ft,(Re).  where  fa('  •  •)  and  ^fe(•  •  •)  have  to 
be  somehow  postulated  and  their  collective  peculiar  consequences  compared  with  available 
data.  The  f|,(-  •  •)  and  Ft,{-  •  )  proposed  here  are  such  that 

b^-2a^=Re> 

and 

6'-H2a^  =  Re-’, 
viz. 

a=0.63(Re-'-  Rfr')"” 
and 

b=0.71(Re-'-^Re>)“^ 


3  APPLICATION  TO  GRID-TURBULENCE 

The  formulated  theory  was  formally  applied  to  the  grid-turbulence  results  of  Stewart  {loc. 
cit.).  First,  a  test  for  seF  •  milarity  in  his  data  was  performed  by  determining  whether  the 
two-point  transfer  function  for  all  three  of  his  cases  (1.  Re  =  5300;  2.  Re  =  21,200;  3.  Re  = 
42,400;  Re  =  UM/v)  could  collectively  be  reasonably  approximated  by  a  single  function  of 
and  r  (or  combinations  of  and  r)  only,  i.e.  no  terms  in  "r"  or  r  to  a  fractional  power;  recall 
that  for  the  two-point  case,  q  =  0.  The  result  of  that  test  is 

K  I-  mr/Mf 
^  ^  li  +  y(r/.v/)^p’ 

where  the  particular  values  of  (/3,  y)  for  each  of  his  Re  are  specified  in  Figs  2-4.  Thus,  for 
theoretical  analyses  of  grid-decay,  a  viable  approximation  of  two-point  third-order  correlations 
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experimental  cur  ve  1 

fi=-l  48  7=70 


Fig  2.  Triple-velocity  curve,  k(-  ■  ■),  obtained  by  Stewart,  broken  line  is  theoretical  curve  (Re  =  5300). 


experimental  curve2 

P  =  -3.0  7  =  12  0 


Fig.  3.  Triple-velocity  curve,  k(-  ■  ■),  obtained  by  Stewart;  broken  line  is  theoretical  curve  (Re  =  21,200) 


experimental  curves 
/J  =  -7  42  7  =  22  0 


Fig.  4.  Triple-velocity  curve,  k(-  ■  ■),  obtained  by  Stewart,  broken  line  is  theoretical  curve  (Re  =  42,400). 
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is 


k{r,  f)  = 


-2a\r/rif 

[\+b\rhj)'r 


where  (a,  b)  can  be  quickly  determined  from  the  Stewart  data  by  visually  approximating  the 
value  of  dimensionless  length,  (r/rj)  in  the  theory  and  {r/M)  in  the  data,-  at  which  k{r/r]) 
reaches  its  maximum  value,  viz.  at  {rlr])  =  \fhb,  and  similarly  the  corresponding  maximum 
value  of  k{r/r]),  viz.  = /c(\/3/6)  =  0.65(a/Z))\  For  the  three  available  cases  these  values 
are  accordingly 


providing  the  values  of 


'0.43,  Re  =  5300 
(u/b)  =  -  0.41.  Re  =  21,200. 
,0.39.  Re  =  42,400 


Note  that  although  b,  through  the  relation.  A:(\/3/6)  explicitly  defines  that  point  along 

the  abscissa  where  k{-  ■  •)  peaks,  through  the  “normalized”  expression, 


K(brlT])  =  ^ibla)-k{brlr}) 


{r/mf 

n  +  (r/r,,)Y 


=  v/b, 


b  also  defines  how  the  length  scale,  rj,  which  is  inherently  peculiar  to  /(•  •  •)  only,  is  itself 
“scaled”  by  the  Re-dependent  b  to  produce  the  length  scale,  rj,,  peculiar  to  k{-  •  •).  Also  note 
that  the  series  expansion, 

K{br/J])  =  -{br/rif[\  -  {br/jjfY  +  ' '  ' 

is  not  bilinear  in /(r/tj)  and  is  in  fact  quadratic  in  f{rlr])  only  for  b  -  0.7  which  as  suggested  by 
the  above  data  is  true  only  for  extremely  low  Rc.  Bilinearity  in  /(•  •  ■)  of  the  leading  term  in  an 
infinite  series  expansion  of  k{-  ■  -)  has  been  (rather  authoritatively)  claimed  to  the  author  by 
Kraichnan  [9].  Note  further  the  virtual  Re-independence  of  {a/b),  a  feature  already  suggested 
by  the  scale-independence  of  C^  (cf.  equation  (3));  this  feature  and  its  implications  will  be 
addressed  later.  The  appropriate  values  of  >•  and  z  in  the  postulated  Fa{-  ■  •)  and  F,,(-  •  •)  are 
then  .v  =  0.17  and  z=0.25,  so  that  for  any  Re  (for  grid-turbulence  only)  a  =  F,(Re)  = 
0.63(Re'’'‘'- Re"'’)'”'  and  6  =  F^fRe)  =  0.71(Re"’' -b  Re"  ”)'”.  .Actually,  the  specific  values 
of  y  and  z  for  each  case  are; 

Re  =  5300;  y  =  0.17.  z=0.25 

Re  =  21,200;  y  =  0. 17,  z  =  0.25 

Re  =  42,000;  y  =  0.18.  z=0.27; 

these  values  suggest  that  z  =  1.5y.  Values  of  the  flow  constants  as  defined  by  Fj(Re)  and 
Ffc(Re).  as  well  as  the  ratio  of  a  over  b,  for  representative  values  of  Re  are  tabulated  in  Table 
1.  The  therein  specified  {a/b)  are  "counter-intuitive"  in  the  sense  that  knowing  a  priori  there 
are  two  modes  manifesting  in  the  flow,  the  natural  inclination  is  to  suspect  some  sort  of 
"equi-partition"  among  them  to  the  extent  that  ideally  a-b  ox  equivalently  {alb)~\. 


Table  1 


Rc 

lO" 

lO' 

10" 

lo’ 

-  j 

10’' 

a 

1  .18 

1  72 

2.14 

264 

3  25 

b 

.154 

4  61 

6  01 

7  87 

10  33 

alb 

0  39 

0  37 

0  36 

0  34 

031 
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Assuming  through  that  a  characteristic  value  of  {alb)  is  roughly  =  0.4,  the  corresponding  value 
of  C*  is  then 

C*  = -8(a/f))’  [  T^do)  = -2.t(o/6)’ = -0.  13,t. 

Jo  (1  +  w-)- 


4.  CONCLUDING  REMARKS 


In  concluding,  it  remains  to  examine  some  of  the  fundamental  implications  of  the  formulated 
theory.  First,  the  complete  solution  of  the  r-integrated  Karman-Howarth  equation,  equation 
(2).  for  grid-turbulence  is  now  possible.  With  the  empirically  determined  value  of  C*  =  -0.41 
in  hand,  the  solution  for  o(r)  =  y,/~"^  and  r}(r)  =  A(t)  =  y;/" '  is  determined  through 
(y,/yi)C,,  = -2Cji  -  16C/ Re"‘.  Re  =  aj//v.  Acceptable  values  of  and  Q  are  fourd  by 
assuming /(r/jj)  =  exp{-»-V27]'},  yielding  -  -Cf  =  \/7r/2  =  5/4;  accordingly. 

5  Re  \ 
vT  Re  +  80/ 


where  y2  =  (10v)"^  As  to  the  constant.  Q,  and  its  virtual  insensitivity  to  changes  in  Re.  the 
physical  explanation  of  this  behavior  appears  to  be  that  even  though  the  “scales"  of  the  largest 
turbulence  eddies  initially  created  by  a  turbulence-generating  source  are  roughly  of  the  same 
order  to  magnitude  of  the  physical  dimensionst  of  said  source  (sort  of  like  a  “spatial- 
frequency”  resonance),  the  particular  mechanism  through  which  these  eddies  are  “scaled 
down”  and  directed  to  the  sink  provided  by  viscous  dissipation  is  defined  not  by  the  source, 
sink,  or  fluid,  but  rather  by  the  mechanisms  of  turbulence  itself;  i.e.,  the  turbulent  fluid,  to  a 
large  extent,  exhibits  what  could  be  described  as  “energy-cascade  autonomy”  or  “self- 
governance”.  Keep  in  mind  that  the  fluid  has  a  strong  inherent  tendency  to  return  to  the  state 
in  which  it  existed  before  the  turbulence  was  created,  and  that  therefore  the  artificially 
produced  eddies,  now  in  a  state  of  excitement,  are  eager  to  “unload"  (on  the  smaller  eddies) 
whatever  portion  of  their  energy  they  themselves  cannot  dissipate;  this  “unloading"  process 
apparently  depends  only  on  turbulence  mechanisms.  More  specifically,  when  an  eddy  of 
scale-size  r  transfers  some  of  its  total  energy,  either  by  entrainment  or  deformation,  to  a  variety 
of  eddies  of  smaller  scale-size  q,  dissipates  the  remainder  of  its  energy  (itself  eventually 
disappearing  unless,  of  course,  it  is  fed  energy  from  those  eddies  of  larger  scale-size),  the 
“transfer”  part  of  such  behavior,  and  its  effectiveness,  is  governed  by  turbulence  mechanisms 
only,  and  like  a  Carnot  engine  in  thermodynamics,  is  independent  of  the  "working  fluid”. 

Second,  the  empirically  derived  expression  for  kirlq)  suggests  that 


s{qln,rlr\)== 


-2a^h2{a){rlqf 
[1  +b-h,{a){r/r})'f' 


Note  that  for  small  enough  r,  s{-  ■  •)  =  -2a^h2{a){r/r}f  while  for  larger  r. 


-2a-^/»;(a)(r/T?)-^ 

b^fi\{a){rlq)* 


This  result  establishes  that  only  one  mode  participates  in  the  transfer  of  energy  for  the  small 
scale-structure  of  the  turbulence  while  both  modes  participate  in  the  large  scale-structure.  This 
means  that  even  though  kinematic  self-similarity  is  assumed,  formulated,  and  pursued,  the 
system  (unsurprisingly)  refuses  to  nonetheless  behave  dynamically  self-similar. 

Third,  a  few  remarks  about  self-similarity.  This  rather  idealized  condition  is  tantamount  to 
saying  that  the  two  “scaling-effects"  of  turbulence,  viz.  of  -  ■  ■)  and  rj(-  ■  ■),  are  both  simply 
functions  of  time  alone,  i.e.  a  -  a(/-only)  and  q  -  r/(/-only).  and  that  as  the  turbulence  decays 
all  terms  in  the  dynamical  equation  consequently  maintain  the  same  relative  balance.  While 


*  Actually,  they  arc  somewhat  smaller  owing  to  some  dissipation  of  energy  by  the  generating  source. 
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within  the  scope  of  the  isotropic  theory  the  /-dependence  only  of  turbulence  intensity  may  be 
completely  justifiable,  the  /-dependence  only  of  the  length  scale  is  unfortunately  not.  In  fact 
the  question  of  different  length  scales  for  turbulence,  viz.  the  dissipation  scale,  A.  for  the 
small-scale  structure  and  the  integral  scale.  A,  for  the  large-scale  structure  {non  self  -similarity). 
has  received  considerable  attention  in  the  literature,  and  the  related  issue  of  a  single  universal 
length  scale  has  never  been  fully  resolved.  The  approach  to  a  multiple  length  scales  formulation 
proposed  here  (for  future  investigation)  is  to  instead  regard  r]{-  •  •)  as  a  function  of  both  r  and 
/,  i.e.  r]{-  •  •)~  t){r,  /).  Since  the  length  scale  is  an  even  function  of  r,  a  possible  expression  for 
T/(r,  /)  is  then 

ri{r,  /)=A-t-(A-A)(r/A)-. 

where  A  and  A  are  those  determined  from  the  self-similar  theory,  such  that  for  r«A, 
7j(r, /)  =  A  while  for  r  =  A,  r/(r, /)  =  A.  Note  that  this  formulation  automatically  provides  for 
differing  length  scales  for  those  values  of  A<r  <  A.  Correspondingly.  k{r,  t)—*k{r/t]{r.  /)}  and 

-2a  V/A)-^[l  +  (r/2A)-l 

[\  +  [I  +2{r/2A)^  +  b\rlk)-  -i-  {rnA)Y 

where  A  =  1.25 A. 

Finally,  the  formulated  thoery  is  not  intended  as  a  “closure"  since  the  value  of  had  to  be 
determined  empirically;  a  bonafide  closure  to  the  turbulence  problem  will  for  self-similar  decay 
theoretically  predict  the  value  of  C*. 
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It  is  theoretically  established  that  energy  transfer  in  nonhomogeneous  turbulence  occurs  in  two 
distinct  scale-independent  “modes.”  The  physical  significance  of  these  modes  is  examined. 
Further,  evidence  is  presented  that  suggests  the  existence  of  a  conservation  principle  that 
governs  the  transfer  of  energy  between  turbulence  eddies  of  differing  size.  A  possible  algebraic 
form  of  such  a  principle  is  formulated. 


Recent  work'”’  has  produced  evidence  suggesting  that 
there  exists  in  self-similar  isotropic  turbulence  a  “law”  that 
governs  the  transfer  of  energy  from  an  eddy  of  one  scale  size 
to  an  eddy  of  smaller  scale  size.  This  law  is  independent  of 
momentum  conservation  (Navier-Stokes  equation)  and  ful¬ 
ly  consistent  with  the  “traditional  wisdom”  that  a  turbulent 
fluid  does  not  transfer  energy  in  an  arbitrary  fashion  but 
indeed  does  so  in  some  well-prescribed  manner.  Recall  that 
in  turbulence,  because  it  is  a  nonlinear  phenomenon,  eddies 
of  differing  (but  nonetheless  comparable)  size  do  indeed 
“communicate”  with  one  another,  and  their  respective  beha¬ 
viors  can  accordingly  not  be  treated  independently.  The  im¬ 
portance  of  such  a  law  to  the  field  of  turbulence,  in  general,  is 
twofold.  Rrst,  it  is  widely  held*  that  the  transfer  of  energy  in 
a  system  obeying  nonlinear  equations  is  a  dynamical  process 
(completely)  dependent  on  those  equations;  yet  the  results 
of  Refs.  1-3  suggest  that,  at  least  in  the  self-similar  isotropic 
case,  there  are  some  aspects  of  energy  transfer  that  are  strict¬ 
ly  independent  of  said  equations  and  in  fact  answer  to  some 
other  (as  yet  unformulated)  calling.  Second,  once  formulat¬ 
ed,  such  a  law  will  provide  a  fundamental  method  for  speci¬ 
fying  how  turbulence  energy  cascades,  and  thus  supply  a 
unified  “closure”  for  the  analysis  of  turbulence  decay. 
Specifically,  the  available  evidence  suggests  that  when  an 
eddy  of  scale  size  r  transfers  energy  to  an  eddy  of  scale  size  q, 
there  is  a  universal  feature  of  said  transfer,  which  is  fully 
independent  of  both  q  and  r,  and  dependent  only  on  the  ratio, 
a  =  q/r,  of  scale  sizes.  This  feature  manifests  for  the  ideal 
self-similar  isotropic  case  in  two  distinct  “modes”  of  energy 
transfer,  denoted  here  as  A,  (a)  andA„  (a),  which  are  joint¬ 
ly  analogous  to  iht  entrainment  /deformation  effects  typical¬ 
ly  observed’  between  eddies  in  traditional  experimental 
studies  of  turbulence,  e.g.,  grid  turbulence,  jet  turbulence, 
etc.  Accordingly,  it  is  not  inappropriate  to  envision  such 
modal  behavior  of  turbulence  energy  transfer  in  some  sense 
analogous  to  the  different  modes  of  classical  heat  transfer, 
save  that  in  turbulence  there  are  only  two  modes  whereas  in 
heat  transfer  there  are  three.  This  point  is  made  here  in  order 
to  emphasize  the  complexity  of  the  overall  turbulence  ener¬ 
gy-transfer  mechanism,  and  perhaps  suggest  the  need  to  ini¬ 
tially  investigate  the  said  mechanism  mode  by  mode  only 
rather  than  in  both  modes  simultaneously.  Kraichnan  and 
Spiegel,*  similarly,  have  already  made  the  analogy  between 
the  transfer  of  kinetic  energy  in  turbulence  and  the  radiative 
transfer  of  electromagnetic  energy  in  an  inhomogeneous  me¬ 
dium  that  both  reflects  and  absorbs.  The  existence  of  such  a 


law  is  predicated  on  the  fact  that  during  self-similar  decay 
the  “scaling  down”  process  of  a  turbulent  eddy  is  invariably 
anal)  tically  modeled  by  introducing  a  length  scale  q  into  the 
decay  formulation  such  that  the  dimensionless  variable  q/q 
in  effect  kinematically  captures  and/or  reflects  the  evolving 
state  of  an  eddy  from  its  inception  (say,  at  the  grid  or  jet)  to 
its  demise  at  that  point  where  the  fluid,  as  a  result  of  viscous 
dissipation,  returns  to  its  original  unagitated  condition. 
(The  distance  of  this  point  from  the  grid  or  jet  is  usually 
referred  to  as  the  “decay  length”  L,  and  L  is  converted  into  a 
“decay  time’^Tby  scaling  out  a  characteristic  flow  velocity, 
viz.,  T=iL/U,  where  centerline  mean-flow  velocity.) 
Note  that  in  any  such  dimensionless  formulation 

a  =  q/rs(q/q)/{r/q), 

regardless  of  whether  q- dissipation  scale  or  q~  integral 
scale,  etc.  The  natural  scale  independence  of  a  automatically 
translates  for  turbulence  into  Re  independence,  and  to  the 
turbulence  practitioner  any  feature  of  fluid  turbulence  that  is 
independent  of  Re  immediately  signals  properties  of  the  tur¬ 
bulence  that  are  unique  and  peculiar  to  turbulence  mecha¬ 
nisms  themselves  and  independent  of  the  physical  character 
of  the  fluid.  In  a  practical  scenario  where  the  turbulence 
under  consideration  is  usually  in  the  steady-state  condition, 
these  unique  and  peculiar  properties  are  also  independent  of 
reference-frame  origin  (homogeneous),  even  though  the 
turbulence  under  investigation  may  itself  be  reference-frame 
origin  dependent  (nonhomogeneous).  Such  "homogeneity 
within  nonhomogeneity”  merits  a  special  place  in  turbulence 
phenomenology  and  deserves  a  careful  investigation  of  its 
ultimate  ramifications. 

.Although  formulated  here  for  turbulence  in  general  for 
the  first  time,  the  initial  motivation  for  such  an  endeavor  was 
furnished  by  Deissler  when  he  concluded’  that  even  for  a 
general  nonhomogeneous  turbulence,  cenain  terms  in  the 
two-point  stochastic-averaged  Navier-Stokes  equation  can 
be  clearly  interpreted  as  “transfer”  terms.  He  accomplished 
this  by  introducing  a  third  point,  defined  as 

x„  =  X  +  n(x' -  X),  0<n<l, 

into  the  two-point  geometry  of  the  flow,  and  cleverl)  deduc¬ 
ing  that  the  appropriately  formulated  “turbulence  self-inter¬ 
action”  terms  produce  zero  total  contribution  to  the  decay  of 
{u,u'),u,-u,(x,t),  u'j-u^ix',t),x'  =x  -t-  r.  The  impact  of 
this  result  (duly  pointed  out  by  Deissler)  is  that  turbulence 
self-interaction  plays  the  same  role  in  nonhomogeneous  tur- 
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*  bulence  that  it  does  in  homogeneous  turbulence,  viz.,  to 
transfer  energy  from  one  eddy  to  another  (smaller)  eddy 
without  changing  the  total  amount  of  mechanical  energy 
present  in  the  flow.  Furthermore,  he  makes  clear  that  the 
interpretation  of  self-interaction  as  “energy  transfer"  is 
accordingly  conceptually  independent  of  the  condition  of 
homogeneity,  recall*  that  the  condition  of  homogeneity 
is  typically  invoked  in  making  the  energy -transfer  inter¬ 
pretation.  More  fundamentally,  though,  Deissler’s  result 
suggests  that  the  two  types  of  turbulence,  through  the  self¬ 
interaction  mechanism,  are  somehow  "indistinguishable.” 
In  this  Brief  Communication  a  correspondingly  similar 
scheme  is  employed  to  show  that  for  nonhomogeneous  tur¬ 
bulence  the  algebraic  form  of  the  self-interaction  tensor, 
Qiji(  ')  ~  when  properly  formulated,  can  also  be 

characterized  in  terms  of  two  invariant  “modes”  of  energy 
transfer,  which  are  identical  for  both  homogeneous  and  non¬ 
homogeneous  turbulence  alike.  The  impact  of  this  new  result 
is  that  it  strengthens  the  unity  perceived  in  the  Deissler 
scheme  and  suggests  the  existence  of  an  energy -transfer  law 
for  all  types  of  turbulence;  the  latter  implies  that  the  premier 
results  of  Refs.  1-3  are  indeed  a  genuinely  unique  feature  of 
turbulence,  and  not  an  artifact  of  the  attendant  simplified 
theory. 

For  simplicity  consider  only  then  component  of  a  three- 
dimensional  turbulence  velocity  field,  and  further  restrict 
the  analysis  to  the  case  where  x  =  (Jt,0,0),  where  .x  —  down¬ 
stream  distance  from  the  grid  or  jet.  In  its  simplest  form, 
then,  the  energy-transfer  mechanism  of  turbulence  is  most 
effectively  modeled  through  the  three-point  correlation 

the  /  dependence  of  C(  •  •  • )  is  suppressed  by  restricting  the 
analysis  to  steady  state  turbulence  only.  Note  that  by  intro¬ 
ducing  the  separation  variables,  g  =  Xj  —  x„  r  =  Xj  —  x„ 
and  the  spatial  location  variable,  3c  =  (x,  -t-  x,  +  Xj)/3, 
2(  ■  •  • )  immediately  converts  into  S{x,q,r).  Through  its  de¬ 
pendence  on  q  and  r,  5(  ■  •  • )  analytically  describes  how  a 
turbulence  eddy  of  size  q  "interacts  with"  (transfers  energy 
to)  an  eddy  of  size  r  or  vice  versa;  i.e..  Si  "  )  reveals  the 
strength  of  the  nonlinearity  of  turbulence  by  quantifying 
how  turbulence  eddies  of  differing  size  communicate  with 
one  another.  In  studies  of  turbulence,  it  is  customary  to  in¬ 
troduce  a  velocity  scale  (typically  a  rms  value)  defined  in 
practice  as 


(7{jc)  =  (-|:j  u'ix,t)dt^  . 


This  scale  permits  0(  ■  •  - )  to  be  naturally  w  ritten  as 
Q(x„x:,x,)  =  a(x,)a(.t,)a(.r,) 

X<c(x„r)t(x,,r)5(xj,/)). 

where  aix)^ix,t)  =  u(x,/),  or  in  terms  of  ix.q.r)  as 
Six,q,r)  —  aix  -  q/l  —  r/3)a(x  -f  Iq/l  -  r/3) 

X<7(x  -  q/3  +  lr/l)six,q,r), 

where  si  ")  =  {^ix^,t)^ixi,t)^ixy,t))  is  a  normalized 
three-point  correlation  and  x,  =  x  -  q/l  -  r/3,  .x,  =  3c 
+  Iq/l  —  r/3,  and  Xj  =  x  -  q/l  -i-  2r/3.  Introducing  the 


correlation  length  scale,  Tj  —  qix),  results  in 

Six.q,r)^aix  ~  q/l  -  r/3)<7(3c  +  2q/l  -  r/l) 

Xc7(x  -  q/l  +  2r/3)r[^/7;(.c),  r/7/(x)), 

and  because  of  the  explicit  presence  of  q  and  r  in  a(  ),  the 
spatial  evolution  of  Si  )  is  not  of  a  self-similar  nature,  seK- 
similarity  follows  only  when  the  length  scale  that  "scales” 
oi  )  is  “large,"  in  comparison  to  rjix).'  Note  that  the 
form  of  velocity  scaling  employed  here  is  quite  different  from 
that  reported  m  Refs.  9  and  10,  in  the  sense  that  for  q^O, 
viz.,  0(x,,x,,x,)  =  <n-(x,)u(x, -f  r)>,  S(x,0,r) 

=a  -(X  -  r/l)aix  -i-  2r/3)^{r/77(x)}  and  not 
a'ix,)aix,  -I-  r)A.(r/7;(x,)),  as  implied  therein.  The  form 
employed  here  duly  reflects  the  maximum  possible  amount 
of  symmetry  (w.r.t.  x„.c,,x,)  in  the  definition  of  x  while 
simultaneously  retaining  the  classical  definitions  of  q  and  r. 
The  rather  “unsymmetric”  velocity  scaling  of  Refs.  9  and  10 
is  perhaps  why  their  triple  correlation  curves  do  not  pass 
through  the  origin  of  coordinates,  and  are  in  fact  all  shifted 
to  the  left  to  the  origin,  thus  suggesting  a  unilaterally  positive 
value  for  A(0).  A  significant  nonzero  value  of  A(0)  means 
that  the  probability  density  function  (pdO  of  u(x,r)  is  mar¬ 
kedly  skewed  about  u  =  0,  and  is  accordingly  inconsistent 
with  the  known”  approximately  normal  distribution  of  tur¬ 
bulence  velocity  at  any  one  point  in  the  flow.  The  scaling  of 
Refs.  9  and  10  follows  from  that  employed  herein  by  impos¬ 
ing  the  transformations  x-x,  +  r/l  for  the  rms  scale  and 
x-x,  for  the  length  scale. 

It  is  a  simple  property  of  Q(  -)  that  Q(x„X2,x,) 
=  Q(  ’fi.-*>-’:2)  =  !2(X2rX|,Xj),  or  equivalently  that  Six,q,r) 
=Six,r,q)  sSix.  -  q,r,  —  q).  Since  the  product  Ic7(x  —  q/ 
3  -  r/l)(rix  -f  2q/l  -  r/l)aix  -  q/l  -f  2r/3)  ]  automati¬ 
cally  satisfies  this  symmetry,  it  follows  thats(  -  •  ■ )  itself  must 
also  satisfy  the  same  symmetry,  i.e.,  six,q,r)=six,r,q) 
=six,  ~  q,r  —  q).  If  si")  is  thus  expanded  in  a  Taylor’s 
scries  in  q  and  r,  the  result  is 

six.q,r)  =:5(x,0,0)  -f  a,ix)q  -f  b,ix)r  -f  a,(.x)9* 
-i-c-ix)qr  +  b2ix)r  +  ■■■  , 


where  the  a„ix),  b„ix),  c„(x),  n  =  1,2,3... ,  represent  the 
appropriately  evaluated  partial  derivatives  of  six,q,r).  It 
can  be  readily  shown'*  ”  that  because  of  the  required  q,r 
argument  symmetry  in  six.q,r),  aside  from  s(.x,0,0),  which 
automatically  satisfies  the  symmetry  the  only  algebraic 
forms  that  can  appear  in  the  Taylor's  expansion  are  the  ho¬ 
mogeneous  polynomials,  Fix,q,r)  =  a-Jx)iq'  -  qr  -r  r) 
and  G(x.g,r)  =  ayix)i2q^  —  Iq'r—  Iqr  -f  2r),  and  those 
higher-order  homogeneous  polynomials,  which  can  be 
formed  from  Fi")  and  G(--),  viz.,  F'i"). 
Fi")Gi"),G'i"),  etc.  No  other  fundamental  algebra¬ 
ic  form  different  from  Fi  ")  and  G( ■  •  ■ )  satisfies  the  re¬ 
quired  symmetry  and  can  therefore  not  appear  in  the  ex¬ 
pansion.  Specifically,  the  9,r dependence  of  six,q,r)  can  thus 
be  fully  characterized  by  Fi  ")  and  Gi  ")  only,  i.e., 
six.q.r)  ~s{Fix,q,r),  Gix,q,r)}.  From  purely  dimensional 
grounds,  then,  it  follows  that  a,(.x)  =  b'/ij^x)  and  flj(x) 
a'/q'ix),  where  a  and  b  are  dimensionless  numbers, 
hence 
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s{  -  }-s{b‘(q-  -qr+r)/v'(x), 

la^iq^  -  \.5q'r  -  l.S^r  +  r^)/q^(.x)}. 

The  physical  significance  of  F(  ')  and  G(  -  )  is  estab¬ 
lished  by  writing  them  as 

F{x,q,r)  =b‘(a'  -  a+  l)r/7/'(x). 

G{x,q,r)  =  2a’(a^  -  l.Sa’  -  1.5a  -t-  \  )r^/qHx), 

and  a  =  q/r,  such  that  a^  -  a -t- 1  sA,  (a)  and 
—  1.5a’-  1.5a -f  l=A„(a)  are  scale-independent  invar¬ 
iant  forms.  Since  rC  }  can  be  fully  characterized  in  terms 
of  the  X  independent  //i(a)  and  A,)  (a),  these  algebraic 
forms  are  then  the  homogeneous  invariant  “modes”  of  non- 
homogeneous  energy  transfer  referred  to  earlier.  The  entire 
nonhomogeneous  transfer  of  energy,  at  <jn>’  spatial  location 
defined  by  x  =  {x,  +  x^  +  x^)/^,  from  an  eddy  of  size 
r  =  x,  -jc„  to  an  eddy  of  size  9  =  x,  -  can  be  construct¬ 
ed  fromA,  (a)  and  A,,  (a),  while  the  related  “spatial  decay” 
of  energy  transfer  leaves  A,  (a)  and  (a)  preserved  and 
unaffected.  The  physical  significance  of  a  and  b  is  thus  that 
they  explicitly  define  “how  much”  of  each  mode  is  present  in 
the  flow;  because  they  are  dimensionless  both  a  and  b  must 
depend  on  Re.  The  mode  interpretation  ofA,  (a)  andA,,  (a) 
is  reinforced  by  introducing  0sa  —  0.5,  whereupon 
b  I  (a)  -/,  W),  h,i  (a)  -/n  (^).  and  noting  that 


for  any  6  (orthogonality).  The  interesting  feature  of  this 
formulation  is  that  depending  on  the  (as  yet  unknown)  ex¬ 
act  algebraic  form  of5(  ■  •  ■ ).  for  “small  enough”  r, 

j(jc,^,r)  =s(x,0,0)  +b^{qr  -qr+  r^)/rjHx) 

~5(3c,0,0)  +  b  ’A,  (a)r/77’(x) 
or 

s(x,q,r)  ~s(3c,0,0)  -f  (a)F/T0{x), 
while  for  “larger”  r 

s{x,q,r}  ~5(x,0,0)  +  b‘h^  (a)r/7]'{x) 

-t-  2a^h,t{a)F/rj-(x)  -f  •••  ; 

this  means  that  for  the  smaller  eddy  sizes  only  one  mode 
participates  in  the  transfer,  while  for  the  larger  eddies  both 
modes  participate.  More  specifically,  even  though  A, (a) 
and  A  „  (a)  are  themselves  “immune”  to  the  effects  of  scale 
size,  there  is  nonetheless  some  mechanism  manifesting  in  the 
flow  which  denies  joint  participating  of  both  modes  in  the 
transfer  ofenergy  for  all  scale  sizes  in  the  flow  This  different 
modal  communication  between  eddies  for  the  large  eddy 
structure  versus  the  small  eddy  structure  has  long  been  sus¬ 
pected  by  turbulence  practitioners,  but  insubstantially  docu¬ 
mented  in  the  literature. 

Since  turbulence  must  transfer  energy  in  a  fashion  that 
preserves  the  scale-independent  forms.  A,  (a)  and  A„(af, 
and  since  Q(  •  •  • )  for  any  turbulence  must  be  algebraically 
expressible  in  terms  of  A  |  (a)  and  A  „  (a),  then  all  universal 
ramifications  of  same  are  embodied  in  these  absolutely  in¬ 
variant  functions,  and  ultimately  in  the  invariant  a.  The 
physical  implication  of  this  resut  is  best  comprehended  by 
recalling  a  simple  experiment  in  physics  which  culminates  in 


a  precise  operational  definition  of  the  mass  of  a  body.  The 
experiment  involves  placing  a  spring  between  two  isolated 
masses,  A  and  B,  pushing  them  together  so  that  the  spring  is 
compressed  ( the  spring  is  not  rigidly  attached  to  either  of  the 
masses),  and  then  releasing  them.  By  repeating  the  experi¬ 
ment  a  number  of  times,  each  time  with  a  spring  of  different 
constant,  it  will  be  found  that  although  the  recoil  velocities 
of  A  and  B  will  be  different  each  time,  the  ratio  /vg  w  ill  be 
identical  for  all  cases,  i.e.,  s  invariant,  regardless  of 
the  amount  of  initial  energy  in  the  system.  The  fact  that  this 
result  provides  a  means  for  unambiguously  defining  the 
mass  of  a  body,  viz.,  mg/m^  — ^nd  that  it  is  consis¬ 
tent  with  the  principle  of  linear  momentum  conservation,  is 
well  established.  Put  simply,  the  invariance  of  Vf/Vg  means 
that  energy  cannot  be  transferred  from  the  potential  energy 
of  the  spring  to  the  kinetic  energy  of  the  masses  in  an  arbi¬ 
trary  fashion,  but  must  indeed  be  transferred  in  a  manner 
consistent  with  linear  momentum  conservation.  In  thermo¬ 
dynamics,  a  correspondingly  simple  experiment  will  estab¬ 
lish  that  in  a  Carnot  engine  the  heat  Q„,  flowing  into  the 
engine  from  a  high  temperature  reservoir  of  fixed  tempera¬ 
ture  ,  is  related  to  the  heat  Qi_ ,  flow  ing  out  to  a  low  tem¬ 
perature  reservoir  of  fixed  temperature  ,  by  the  relation 

Qh/Ql  =  r„/riS  invariant; 

the  invariance  of  this  ratio  is  a  consequence  of  the  principle 
of  conversation  of  entropy,  i.e.,  heat  is  not  transferred  arbi¬ 
trarily  through  the  engine.  In  geometrical  optics  the  fact  that 
the  sines  of  the  angles  of  incidence  and  refraction  stand  in  a 
constant  ratio  to  each  other,  regardless  of  the  magnitude  of 
each,  leads  to  the  formulation  of  Snell's  law.  It  appears, 
therefore,  that  the  invariance  of  a  =  q/r  is  a  manifestation 
of  some  conservation  law  or  principle  which  governs  the 
transfer  of  energy  in  turbulence.  Such  a  conservation  law 
must  describe  the  means  of  (scale-independent)  communi¬ 
cation  between  turbulence  eddies  of  different  size  such  that 
in  the  limit,  eddies  distinguished  by  a  great  disparity  in  size 
should  exert  minimal  influence  on  one  another.  Moreover, 
such  a  principle  is  consistent  with  accepted  wisdom  that  tur¬ 
bulence  does  not  transfer  energy  arbitrarily  but  indeed  does 
so  in  some  well-prescribed  fashion.  If  a  turbulence  eddy  of  a 
particular  size  is  regarded  as  a  dynamical  system  which  re¬ 
ceives  energy  from  those  eddies  of  larger  size,  dissipates 
some  of  this  energy,  and  transfers  the  rest  to  those  eddies  of 
smaller  size,  then  a  possible  such  law  is 

Er{r)/SEg  (r,q)  =  {r/qY-. 

in  this  expression  £r(r)  is  the  total  energy  transferred  by  an 
eddy  ofsizertoall  other  eddies  of  smaller  size,  SEg  (.r.q)  is 
that  portion  of  Ef{r)  that  is  eventually  rece/te</ by  an  eddy 
of  size  y  IS  a  nonadjustable  constant  greater  than  zero,  and 
r>q.  The  unavoidable  x  dependence  of  both  £7-  ( •  ■  ■ )  and 
SEg  [■■  ■)  is  explicitly  deleted  in  the  above  expression  in 
order  to  emphasize  the  location  independence  of  the  ratio 
Er(  -  )/SEg  (  •  • ).  The  total  energy  received  by  an  eddy  of 
size  q  is  thus 

.V 

Eg(q^)  =  ^  SEg{r„,q^x), 

t  =  I 

where  r„^n  —  1,2,3  . ,  .Vare  all  greater  than  q.  The  ultimate 
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'*  consequences  of  such  a  law  (or  others  like  it )  need  to  be  fully 
investigated. 
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APPENDIX  B 


Pressure  and  Velocity  Correlations 
in  a  Plane  Turbulent  Water  Jet 
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